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Abstract 


The  topic  that  was  investigated  was  whether  there  is  a  lift  and  thrust 
benefit  to  flapping  an  underwater,  biologically  inspired  three  dimensional  foil 
near  the  solid  bottom  surface  (so  that  it  was  in  ground  effect).  The  experimental 
method  used  a  dual  canister  device  that  allowed  actuation  in  roll  and  pitch,  with 
force  sensors  attached  to  the  pitch  shaft,  in  order  to  record  the  forces  produced 
by  the  foil  as  it  flapped.  The  dual  canister  was  towed  on  a  carriage  at  a  constant 
speed  in  a  large  tow  tank  that  had  been  configured  to  have  a  long  run  of  constant, 
deep  water  depth,  a  very  short  transition  period,  and  then  another  long  run  of 
constant,  shallow  water  depth.  Due  to  this  configuration,  in  one  run  the  foil  was 
able  to  encounter  freestream  conditions  and  ground  effect  conditions.  The 
results  proved  in  all  sixteen  experimental  cases  for  varying  Strouhal  number  and 
maximum  angle  of  attack  that  the  mean  lift  coefficient  near  the  bottom  was  larger 
than  that  in  the  freestream.  A  potentially  useful  data  point,  for  which  the  mean 
thrust  coefficient  had  a  positive  change  from  flapping  near  bottom  as  compared 
to  the  freestream,  was  found  to  have  a  change  in  maximum  instantaneous  lift 
force  of  14%.  This  would  give  a  large  enough  change  in  signal  strength  that  it 
could  be  used  as  a  parameter  on  a  future  underwater  vehicle  to  control  altitude 
above  the  ground.  The  benefit  to  flapping  in  ground  effect  was  equivalent  to  a 
larger  than  1°  pitch  bias  difference  at  a  zero  mean  lift  coefficient.  Additionally, 
there  was  a  thrust  benefit  seen  to  flapping  in  ground  effect,  but  only  under 
certain  kinematics.  Though  not  as  dramatic  as  the  benefit  in  lift,  there  was  still  an 
8%  difference  in  the  mean  thrust  coefficient  observed  between  flapping  near  the 
bottom  and  flapping  in  the  freestream,  for  the  case  where  the  largest  change  in 


mean  thrust  coefficient  was  observed.  This  could  equate  to  a  large  savings  in 
battery  life,  and  hence  a  longer  endurance  for  a  vehicle  taking  advantage  of  the 
thrust  benefit  seen  by  flapping  in  ground  effect. 

While  this  work  remains  preliminary  in  nature,  it  shows  that  much  more 
useful  work  remains  to  be  done  to  explore  the  benefit  induced  by  flapping  a  foil 
close  to  a  hard  surface  ground. 
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1.  Introduction 


1.1.  Motivation 

Aquatic  animals  that  employ  flapping  fins  for  underwater  propulsion,  such 
as  turtles,  fish,  and  penguins,  are  endowed  with  amazing  abilities  that  engineers 
have  yet  to  even  match  in  our  unmanned  underwater  vehicles  (UUV).  These 
animals  are  able  to  gracefully  swim,  navigate,  avoid  obstacles,  avoid  predators, 
and  catch  prey  in  complex  underwater  environments.  Some  of  these  complex 
underwater  environments  include  areas  of  interest  that  humans  greatly  desire  to 
operate  UUVs  in.  There  are  many  reasons  why  flapping  foil  UUVs  would  need  to 
operate  close  to  the  bottom  surface  of  a  body  of  water  such  as  the  ocean  floor. 
Navigation  would  be  improved  and  the  adverse  effects  of  waves  and  currents 
would  be  minimized.  Military  UUVs  could  take  advantage  of  the  stealth  that  it 
provides,  as  well  as  simply  the  opportunity  to  closely  inspect  long  stretches  of 
pipeline  and  communications  cabling.  In  the  future,  fin  propelled  vehicles  will  be 
able  to  use  real  time  fin  force  data  to  estimate,  and  hence  control,  altitude  near 
the  bottom  surface  (Licht  &  Dahl,  2013).  This  work  is  a  humble  beginning  to 
investigating  the  benefits  and  challenges  that  may  be  observed  by  flapping  in 
ground  effect. 

1.2.  Chapter  Preview 

The  rest  of  this  chapter  includes  background  and  literature  review 
sections,  which  describe  key  terms  related  to  flapping  foils.  Chapter  2  presents 
the  methodology  of  this  work,  including  improvements  that  were  made  to  the 
testing  platform  that  was  used,  and  the  experimental  method.  Chapter  3 
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describes  the  results  that  were  obtained,  presenting  the  data  in  a  variety  of  ways, 
and  discusses  the  findings.  Chapter  4  suggests  further  design  improvements  that 
are  recommended,  and  presents  some  error  analysis.  Chapter  5  proposes  some 
future  work  that  may  be  done  with  the  testing  platform  with  respect  to  ground 
effect.  Chapter  6  gives  the  summary  and  conclusions  from  this  work.  Appendices 
and  a  bibliography  follow  chapter  6. 

1.3.  Background  and  Literature  Review 

1.3.1.  Key  terminology 

1.3. 1.1.  Foil  Dimensions 

The  flapping  foil  that  has  been  used  in  this  work  is  one  of  the  four  that 
were  employed  on  Finnegan  the  RoboTurtle,  described  in  (Licht,  2008).  It  has 
been  made  from  a  titanium  framework  surrounded  by  Shore  80A  polyurethane 
elastomer,  in  a  biologically  inspired  shape  that  emulates  a  turtle  fin.  The  chord 
(c)  is  the  distance  from  leading  edge  to  trailing  edge.  Since  it  is  not  a  rectangular 
planform,  the  chord  that  will  be  used  is  0.1m,  which  is  the  mean  chord.  In  this 
work  the  mean  chord  will  always  be  referred  to  as  simply  "chord".  The  span  (s) 
is  the  distance  from  root  to  tip,  0.4m.  The  cross  section  of  the  fin  is  symmetrical 
about  its  chord,  meaning  that  it  is  not  cambered.  At  its  thickest,  the  foil  is  0.015m 
thick.  The  cross  section  approximates  a  NACA  0012  shape  with  a  rounded 
leading  edge  and  a  sharp  trailing  edge.  Foil  dimensions  are  depicted  below. 
Because  of  the  polyurethane  elastomer  it  is  constructed  from,  and  the  placement 
of  the  titanium  framework,  one-third  of  the  chord  closest  to  the  trailing  edge  is 
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compliant.  This  property  has  not  been  quantified  nor  studied  in  the  current 
work. 


Figure  1.  Flapping  foil  dimensions 


1.3. 1.2.  Non  dimensional  parameters 

It  is  well  known  that  a  foil  oscillating  in  heave  and  pitch  can  produce  a 
mean  forward  thrust  force,  through  the  generation  of  a  reversed  von  Karman 
vortex  street.  This  is  the  principle  by  which  many  swimming  and  flying  animals 
propel  themselves  (Streitlien  &  Triantafyllou,  1997).  The  depiction  below  shows 
a  reversed  vortex  street  generated  behind  a  fish  swimming  by  flapping  its  caudal 
fin  (tail),  as  seen  from  above. 


Figure  2.  Reverse  vortex  street  behind  a  swimming  fish  (Eloy,  2012) 


Where  U  is  the  forward  speed,  A  is  the  peak-to-peak  amplitude  of  the  tail  flap,  and 
b/a  is  the  vortex  spacing  ratio.  This  same  pattern  is  generated  by  the  flapping  of 
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pectoral  fins  on  fish  using  labriform  motion  (Sfakiotakis,  Lane,  &  Davis,  1999), 
and  by  the  fins  on  turtles.  When  averaged  over  one  or  more  periods  of  motion  in 
time,  the  vortex  street  has  the  velocity  profile  of  a  jet,  with  zero  mean  lift  force, 
and  some  non-zero  mean  thrust  force. 

The  first  important  non  dimensional  parameter  is  the  Strouhal  number 
(St),  a  ratio  used  to  characterize  the  vortex  pattern  as  it  relates  the  vehicle 
velocity  to  the  frequency  and  size  of  each  vortex  generated.  In  this  work,  St  has 
been  defined  as: 

f'A  =  /  ■  2r0  7  ■  (po 
U  U 

Equation  1.  Definition  of  Strouhal  number 

Where /is  the  flapping  frequency,  ro.7  will  be  addressed  below,  and  cpo  is  the  roll 
amplitude.  For  a  two  dimensional  foil,  the  amplitude,  A,  is  clear.  However,  for 
the  three  dimensional  case  it  is  less  clear  what  distance  should  be  used. 

Following  (Techet,  2008),  the  amplitude  to  base  St  upon  is  taken  as  the  amplitude 
of  an  arc  length  at  70%  of  the  foil  span,  away  from  the  root.  The  location  was 
selected  to  be  consistent  with  conventional  propeller  notations  and  for  easy 
comparison  with  past  flapping  foil  experiments.  The  value  ro  is  the  distance  from 
the  roll  axis  to  the  root  of  the  fin.  The  radius  of  the  70%  span  location  of  the  foil 
is  denoted  as  ro.7  with  the  roll  axis  as  the  origin: 

?o.7  =  r0  +  0.7s 

Equation  2.  Radius  of  70%  span  location 
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Following  (Polidoro,  2003),  the  location  of  ro.7  is  very  close  to  the  effective  center 
of  hydrodynamic  force  on  the  foil.  For  this  reason  in  the  current  work,  ro.7  is 
assumed  to  be  the  location  of  the  hydrodynamic  center  of  the  foil,  the  location  at 
which  the  lift  forces  and  thrust  forces  generated  by  the  entire  fin  is  taken  to  act 
through.  This  location  is  0.28m  away  from  the  root,  and  0.12m  away  from  the  tip. 
One  final  note  is  that  St  in  Equation  1  is  based  on  the  amplitude  of  ro.7,  and  not 
based  on  the  peak-to-peak  amplitude  of  the  fin  tip. 

The  next  important  non  dimensional  parameter  is  the  heave  to  chord  ratio 
( ho.7/c ).  This  ratio  is  the  heave  amplitude  divided  by  the  chord  length,  and  is 
based  on  the  heave  amplitude  at  the  assumed  hydrodynamic  center. 

Lastly,  height  above  ground  to  chord  ratio  ( H/c )  will  be  used  to 
characterize  when  the  foil  is  in  ground  effect.  For  H/c  of  greater  than  3,  the 
effects  of  being  close  to  the  ground  are  negligible  (Wu,  Shu,  Zhao,  &  Yan,  2014). 

So  by  that  definition,  ground  effect  may  be  expected  for  H/c<3. 

1.3. 1.3.  Foil  Kinematics 

A  few  more  equations  and  figures  must  be  introduced  in  order  to 
characterize  the  foil  kinematics  (Polidoro,  2003).  Roll  motion  is  described  by: 

<p(t)  =  cp  o  sin(<ut) 

Equation  3.  Equation  for  roll  of  the  foil 

Where  co  is  the  flapping  frequency  expressed  in  radians.  Similarly,  pitch  motion 
can  be  described  by: 

0(t)  —  90  cos(wt) 

Equation  4.  Equation  for  pitch  of  the  foil 
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Where  do  is  the  pitch  amplitude.  In  this  work  the  phase  angle  ( ip )  between  roll 
and  pitch  will  always  be  tt/2,  such  that  maximum  pitch  occurs  at  zero  roll,  and 
zero  pitch  occurs  at  maximum  roll.  The  angle  of  attack  encountered  by  the  foil  at 
each  point  along  the  span  has  two  components.  The  first  is  due  to  the  ratio  of 
forward  velocity  and  heave  velocity,  and  the  second  is  due  to  the  instantaneous 
pitch  position  of  the  foil: 


Figure  3.  Angle  of  attack  at  one  span  location  (Polidoro,  2003) 


Therefore  the  angle  of  attack  is  described  by: 

^  *  /W0.7<Po  cos(wt)\  ,  „  r  ^  ,  n 

a(t)  =  —arctan  I - — - I  +  90  cos  (cut)  +  9bias 

Equation  5.  Equation  for  angle  of  attack  of  the  foil 

The  first  term  represents  the  roll  induced  angle  of  attack  and  the  second  two 
terms  represents  the  pitch  induced  angle  of  attack.  For  the  majority  of  motions 
that  produce  thrust,  the  pitch  motion  is  selected  to  reduce  the  maximum  angle  of 
attack  ( ccmax ).  The  angle  of  attack  is  a  function  of  spanwise  location,  but  again  the 
maximum  angle  of  attack  is  a  parameter  that  is  calculated  at  the  assumed 
hydrodynamic  center  of  the  foil,  as  seen  in  the  equation  above. 
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As  a  performance  metric,  the  mean  lift  and  mean  thrust  coefficients  will  be 

used  later  on  in  analysis.  The  mean  lift  coefficient  is  found  by  the  following: 

_  _  2L 
^ L  pU2sc 

Equation  6.  Equation  for  mean  lift  coefficient 

Where  p  is  the  water  density.  Similarly,  the  mean  thrust  coefficient  is: 

2  T 

CT  - - 

pU2sc 

Equation  7.  Equation  for  mean  thrust  coefficient 

For  any  given  geometry  the  St  and  amax,  when  taken  together,  are  the  only  things 
needed  to  completely  describe  the  foil  motion  (Polidoro,  2003). 

1.3.2.  Ground  effect  for  fixed  airfoils 

To  start  the  literature  review  section,  the  first  paper  that  will  be  discussed 
relates  to  a  simple  case  in  aerodynamics  using  a  fixed  airfoil,  moving  close  to  the 
ground.  Ground  effect  is  well  understood  in  aerodynamics  for  fixed  airfoils. 
(Garcia  &  Katz,  2003)  described  trapped  vortices  between  fixed  airfoils  and  the 
ground  surface,  giving  insight  into  their  application  on  racecars.  The  principle  of 
increasing  fluid  dynamic  loads  (lift)  by  creating  strong  vortices  near  solid 
surfaces  existed  in  nature  long  before  attempts  were  made  to  understand  its 
mechanics.  It  is  natural,  therefore,  that  vorticity  is  frequently  used  to  explain  an 
airfoil’s  lift,  which  is  often  called  "bound  vorticity"  (to  separate  it  from  "unbound” 
vortices  found  in  wakes).  Consequently,  augmentation  of  the  fluid  dynamic  loads 
by  adding  unbound  (or  trapped)  vortices  is  a  logical  extension  to  the  "more 
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vorticity,  more  lift”  principle.  The  authors  concluded  that  the  trapped  vortex  is  a 
viable  principle  for  lift  augmentation,  and  immediately  moved  to  apply  that  to 
ground  vehicles,  in  particular  racecars.  Ground  vehicles,  contrary  to  aircraft, 
move  close  to  the  ground,  and  their  incidence  relative  to  the  freestream  falls  into 
the  small  angle  of  attack  (or  zero  incidence)  category.  The  utilization  of  this  type 
of  vortex  lift,  therefore,  becomes  quite  attractive.  For  example,  such  vortices 
when  trapped  beneath  a  moving  automobile  and  the  ground  can  increase  the 
vehicle’s  negative  lift  (downforce).  This  aerodynamic  downforce  is  often 
exploited  by  racecar  designers  in  order  to  increase  tire  adhesion  and  vehicle  high 
speed  performance. 

1.3.3.  Animal  flight  aerodynamics  in  ground  effect 

In  this  next  paper,  the  author  studied  flying  animals,  making  general 
observations  about  their  performance  flying  in  ground  effect.  Again,  fixed  wings 
were  used  in  this  study.  (Rayner,  1991)  developed  a  complete  theory  for  a  fixed 
wing  in  ground  effect,  based  on  a  steady  state  lifting-line  wing  model,  and  uses  it 
to  show  how  ground  effect  is  likely  to  affect  the  flight  performance  of  an  animal 
close  above  a  surface.  The  most  important  findings  were  that  flight  in  ground 
effect  above  a  flat,  smooth  surface  may  give  an  animal  considerable  performance 
advantages,  including  a  reduction  in  cost  of  transport  of  up  to  15%,  and  a 
reduction  in  mechanical  flight  power  of  as  much  as  35%,  compared  with  values 
for  flight  out  of  ground  effect.  The  author  additionally  concluded  that  slow  flight 
performance  in  ground  effect  is  very  poor,  owing  to  the  horizontal  air  velocity 
induced  around  the  wing  in  the  presence  of  the  ground.  Ground  effect  was 
defined  in  that  work  as  the  situation  when  an  animal  or  an  aircraft  flies  close 
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above  a  plane  surface;  the  aerodynamic  properties  of  its  wings  are  altered  as  a 
result  of  the  interaction  of  the  vortices  on  the  wing  and  the  wake  with  the 
surface.  Ground  effect  was  noted  to  have  a  significant  effect  on  flight 
performance  during  takeoff  and  landing.  The  dominant  effect  is  a  reduction  in 
the  induced  (vortex)  drag,  and  hence  a  saving  in  the  thrust  required  for  level 
flight.  In  the  current  study  it  would  be  interesting  to  find  a  similar  thrust  benefit. 

1.3.4.  Ground  effect  on  an  infinite  foil 

Moving  away  from  both  fixed  wings  and  aerospace  applications,  the  next 
paper  conducted  preliminary  experimental  studies  in  force  production  by 
heaving  and  pitching  oscillating  (two  dimensional)  foils  in  ground  effect  in  a 
small  water  tow  tank  (Licht  &  Dahl,  2013).  They  showed  that  the  mean  lift  for  a 
typical  set  of  thrust  generating  foil  kinematics  when  operating  within  two  chord 
lengths  of  the  bottom  is  increased  by  an  amount  consistent  with  approximately  a 
1  degree  positive  (away  from  wall)  bias  in  the  foil  pitch  angle.  Additionally,  they 
found  that  flapping  a  foil  near  a  solid  boundary  generates  significant  (18%) 
variation  in  peak  magnitude  of  downstroke  vs.  upstroke  instantaneous  lift.  In 
this  work,  there  is  a  dual  canister  mechanism  that  enables  the  foil  motion.  This 
also  closely  simulates  a  vehicle  body.  In  Licht  and  Dahl’s  work  only  the  foil  was  in 
the  flow  and  the  driving  mechanism  was  outside.  Additionally,  their  foil  span  was 
vertical  and  flapped  near  a  wall,  but  in  this  work  the  foil  is  flapped  near  a  ground. 
This  is  bound  to  create  differences  in  the  magnitude  of  the  effect  that  is  observed 
in  the  two  data  sets. 

1.3.5.  Fluid  dynamics  in  ground  effect 
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The  next  work  numerically  studied  insect  wing  flapping  in  air.  (Wu,  Shu, 
Zhao,  &  Yan,  2014)  investigated  ground  effect  on  a  flapping  insect  wing  in 
forward  flight,  through  simulation  using  an  Immersed  Boundary-Lattice 
Boltzmann  Method  (IB-LBM).  A  NACA  0012  airfoil,  which  models  the  insect  wing 
cross-section,  was  considered.  The  airfoil  executed  a  combined  motion  of 
harmonic  heave  and  pitch  rotation.  The  Reynolds  number  and  the  amplitude  of 
motion  were  fixed,  while  the  height  above  ground  and  frequency  of  oscillation 
were  examined.  They  found  that  the  flow  patterns  shed  from  the  foil  were 
altered  due  to  the  ground  effect.  The  following  two  figures  show  the 
instantaneous  vorticity  contours  under  different  Strouhal  numbers.  Looking  at 
flow  patterns  alone,  the  authors  concluded  that  there  was  negligible  difference  at 
an  H/c>3,  which  is  used  as  the  basis  for  freestream  conditions.  In  these  figures, 
the  solid  line  represents  counterclockwise  vorticity  and  the  dashed  line 
represents  clockwise  vorticity. 


Figure  4.  Instantaneous  vorticity  contours  for  H/c=3  (freestream) 
(Wu,  Shu,  Zhao,  &  Yan,  2014) 
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Figure  5.  Instantaneous  vorticity  contours  for  H/c=l  (near  ground) 
(Wu,  Shu,  Zhao,  &  Yan,  2014) 


Note  that  the  vortex  interaction  with  the  ground  is  very  significant  at  H/c=l, 
which  changes  the  vortex  structures  in  the  wake.  At  low  St  (St=0.1  and  0.2),  the 
size  of  the  vortex  shed  from  the  flapping  foil  is  shortened  due  to  the  presence  of 
ground.  However,  the  strength  of  vortices  is  increased,  which  may  be  attributed 
to  vortex  interaction  with  the  ground.  As  St  keeps  increasing  (St=0.3  and  0.4), 
the  vortex  interaction  becomes  stronger  and  stronger,  gradually  affecting  vortex 
shedding  of  the  foil.  When  St=0.5,  it  is  most  evident  that  the  vortices  have  been 
compressed  to  an  oblate  shape.  Consequently,  there  is  a  notable  angle  between 
the  ground  and  the  center  line  of  the  vortex  street.  The  authors  relate  this  angle 
to  changes  in  the  overall  mean  lift  vector  direction,  resulting  in  an  increased 
mean  lift  coefficient  seen  in  ground  effect.  The  authors  also  noted  that  the  mean 
drag  coefficient  was  increased  at  low  frequency  (low  St)  and  decreased  at  high 
frequency  (high  St),  i.e.,  the  mean  thrust  coefficient  would  decrease  for  low  St  and 
increase  at  high  St. 

1.3.6.  Swimming  near  the  substrate 

Having  just  looked  at  a  numerical  simulation  which  investigated  ground 
effect,  the  next  report  to  be  introduced  is  the  most  relevant  to  the  present  work. 
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In  (Blevins  &  Lauder,  2013)  a  simple,  stingray-inspired  physical  model  was  used 
to  give  insight  into  ground  effects  on  undulatory  swimmers,  contrasting  the  self- 
propelled  swimming  speed,  power  requirements,  and  hydrodynamics  of  fins 
swimming  with  fixed  kinematics  near  and  far  from  a  solid  boundary.  The  authors 
experimentally  determined  that  contrary  to  findings  for  gliding  birds  and  other 
fixed-wing  fliers,  ground  effect  does  not  necessarily  enhance  the  performance  of 
undulating  fins.  Under  most  kinematic  conditions,  fins  do  not  swim  faster  in 
ground  effect,  power  requirements  increase,  and  the  cost  of  transport  can 
increase  by  up  to  10%.  The  influence  of  ground  effect  varies  with  kinematics, 
suggesting  that  benthic  fish  might  modulate  their  swimming  behavior  to 
minimize  locomotor  penalties  and  incur  benefits  from  swimming  near  a 
substrate.  In  this  work  a  recirculating  tank  was  used  and  the  flow  speed  was 
altered  to  match  the  thrust  produced  by  the  moving  fin.  When  the  fin  maintains 
an  equilibrium  position,  thrust  and  drag  are  balanced  during  each  cycle  of 
motion.  Although  the  experimental  setup  was  different,  the  results  from  this 
study  will  prove  to  be  valuable  in  the  current  work. 

1.3.7.  Testing  platform 

The  current  effort  extends  and  builds  on  (Rauworth,  2014)  wherein  a 
system  capable  of  testing  submerged  underwater  flapping  foils  in  a  tow  tank 
while  recording  force  and  position  data  was  created.  The  existing  testing 
platform  consisted  of:  a  dual  canister  system  to  enable  roll  and  pitch  of  a  foil,  an 
aluminum  carriage  attachment  assembly,  an  instrument  chassis,  a  power  supply, 
a  laser  distance  measurement  (LDM)  device  on  the  tow  tank  carriage,  and  a 
computer  on  the  shore  to  both  control  the  dual  canister  system  and  to  record  the 
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data  from  the  force  sensors.  It  is  this  testing  platform,  as  well  as  the 
experimental  techniques  and  processing  codes  that  the  author  developed,  that 
has  been  used  in  the  current  work.  The  following  two  figures  show  pictorial 
representations  of  the  testing  platform: 


Carriage 
Platform  With  ■ 
Components 


ta5 


Computer 
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Tank 


-  Carriage 


_  Fiber-optic 
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Figure  6.  Plan  view  of  testing  platform  (Rauworth,  2014) 


DAQ 

Board 


Figure  7.  Profile  view  of  testing  platform  (Rauworth,  2014) 
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In  addition  to  setting  all  of  the  above  listed  components  up  and  getting  them  to 
work  together,  the  author  also  compared  his  newly  generated  data  to  that  of  past 
testing  conducted.  For  a  comparison  of  mean  thrust  coefficient  contour  plots, 
Rauworth  used  (Techet,  2008).  In  that  work,  a  linearly  tapered  three 
dimensional  foil  was  used  with  the  same  dual  canister  apparatus,  a  recirculating 
tank,  and  a  six-axis,  waterproof  strain  gauge  sensor.  The  sensor  measured  the 
three  force  components  and  the  three  moment  components  on  the  foil.  Measured 
forces,  mean  force  coefficients,  and  hydrodynamic  efficiency  data  were 
presented.  Mean  thrust  coefficient  contours  from  (Techet,  2008)  and  (Rauworth, 
2014)  closely  matched  in  shape  and  trend,  but  the  values  were  different  (due  to 
different  foils  and  sensing  methods  used).  For  a  comparison  of  time  sequenced 
lift  and  thrust  forces,  Rauworth  used  (Polidoro,  2003).  The  two  displayed  similar 
trends  and  waveforms  when  compared.  Rauworth  found  that  phase  averaged 
data  adheres  to  expected  theoretical  results  and  a  representative  figure  is  shown 
below,  for  St=0.5  and  amax=30°.  It  is  included  to  summarize  the  author’s  findings 
and  as  an  example  to  compare  the  results  of  the  current  work  with  in  a  later 
section: 
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30 


Phase(degrees) 

Figure  8.  Phase  averaged  lift  and  thrust  plot  example  (Rauworth,  2014) 


Finally,  in  chapter  5  of  (Rauworth,  2014),  he  suggested  design  improvements  that 
should  be  made  to  the  testing  platform,  and  virtually  every  component  has  been 
improved  upon  in  the  current  work. 
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2.  Methodology 

2.1.  Testing  Platform  Improvements 

The  following  figure  shows  the  body  frame  coordinate  system  of  the  dual 
canister.  Notice  that  the  origin  of  the  coordinate  system  is  established  at  the 
intersection  of  the  roll  axis  and  the  pitch  shaft. 


Figure  9.  Body  frame  coordinate  system 


2.1.1.  Install  Force  Sensors  with  Larger  Separation 

The  two  Kistler  type  9602  three-axis  force  sensors  are  attached  via 
bearings  to  the  pitch  shaft.  The  lift  and  thrust  forces  on  the  foil  are  transmitted  to 
the  force  sensors,  which  send  the  recorded  forces  (sampled  at  a  rate  of  200Hz)  to 
the  shore  based  computer.  The  procedure  to  use  the  recorded  forces  in  all  three- 
axes  to  find  the  lift  and  thrust  on  the  foil  is  described  later  in  the  Data  Processing 
section.  The  force  sensors  were  previously  installed  at  a  y-distance  separation  of 
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0.03m,  due  to  design  constraints.  In  the  current  work,  the  force  sensors  are 
installed  at  a  y-distance  separation  is  0.092m,  a  threefold  increase.  The  desire  to 
increase  the  separation  between  the  force  sensors  was  driven  by  the  need  to 
improve  the  signal  to  noise  ratio. 

2.1.2.  Rigid  Bearings  to  Spherical  Bearings 

The  pitch  shaft  was  previously  supported  by  two  large,  rigid  bearings, 
which  only  allowed  rotation  about  the  pitch  axis.  Acting  as  clamped  connections, 
the  bearings  undoubtedly  imparted  moments  about  the  body  x-  and  z-axes.  A 
representation  follows,  showing  first  the  closely  spaced,  large,  rigid  bearings, 
which  clamped  around  the  pitch  shaft.  Below  that  shows  the  spaced  out,  small, 
spherical  bearings,  which  act  as  pinned  supports  to  the  pitch  shaft. 


Figure  10.  Clamped  vs.  pinned  pitch  shaft  support  bearings 


It  was  imperative  to  eliminate  these  additional  moments,  because  they 
introduced  significant  errors.  The  method  of  obtaining  lift  and  thrust  forces  for 
the  flapping  foil  involves  using  force  data  from  the  force  sensors,  finding  the 
moments  of  these  forces  about  the  x-  and  z-axes,  and  then  dividing  out  a  y- 
distance  to  the  fin  assumed  hydrodynamic  center  (equations  given  in  a  later 
section  in  this  chapter).  The  solution  was  to  replace  the  rigid  bearings  with 
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spherical  bearings,  which  readily  allow  rotation  about  the  pitch  axis,  but  are  not 


constrained  in  rotation  about  the  other  two  axes. 


Figure  11.  Comparison  between  old  (large  metal)  and  new  (small  plastic)  bearings 

The  bearings  chosen  were  igubal™  pillow  block  bearings  from  igus  incorporated. 
The  new  bearings  are  maintenance  free  and  self-lubricating.  The  igubal  bearings 
consist  of  a  housing  with  a  spherical  plastic  insert  which  freely  rotates  in  any 
direction  inside  the  housing.  Although  stiction  may  be  a  problem  with  the  new 
bearings,  they  will  facilitate  a  far  smaller  degree  of  data  corruption  due  to 
extraneous  moments,  compared  to  the  previous  rigid  bearings.  An  additional 
benefit  to  replacing  the  bearings  is  a  significant  space  savings  in  the  pitch 
canister,  allowing  for  easier  installation  of  internal  components. 

2.1.3.  Split  plate  design 

The  previous  method  of  mounting  the  two  large,  rigid  bearings  had  them 
rigidly  mounted  onto  a  solid  aluminum  plate  which  rested  on  top  of  both  force 
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sensors,  which  were  on  top  of  another  solid  aluminum  plate  secured  into  the 
delrin  canister  housing.  The  solid  plate  that  connected  the  force  sensors  also 
potentially  supported  moments,  again  corrupting  the  data.  The  single  plate  was 
replaced  with  two  small  aluminum  plates.  The  new  mounting  system  is  a  solid 
aluminum  plate  secured  into  the  delrin  housing,  with  the  two  force  sensors  above 
that.  On  top  of  force  sensor  1  is  a  small  aluminum  plate  that  supports  one  igubal 
bearing,  the  pitch  motor,  and  the  pitch  homing  bracket.  On  top  of  force  sensor  2 
is  another  small  aluminum  plate  that  supports  the  second  igubal  bearing.  The 
new  mounting  system  eliminates  any  solid  connection  that  bridges  the  gap 
between  the  bearings  and  force  sensors.  Again,  an  additional  improvement  is  the 
space  and  weight  savings  of  the  new  system,  as  well  as  easier  installation. 

2.1.4.  Wiring  Improvements 

2. 1.4.1.  Method  of  connection  of  dual  canister  to  carriage 

Cables  conducting  the  data  from  the  force  sensors,  through  the  dual 
canister,  up  the  carriage  attachment  to  the  carriage,  used  to  attach  to  a  National 
Instruments  block  (a  NI  SCB-68A  shielded  68-pin  connector  block  for  DAQ 
devices)  via  rather  tenuous  four  pin  connectors  that  frequently  bent  or  broke. 
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Figure  12.  Former  force  sensor  data  cable  attachments 


Rather  than  having  two  connections  to  make  each  time  the  system  is  hooked  up 
to  the  carriage,  many  wires  were  eliminated,  and  the  two  cables  were  hard-wired 
directly  into  the  NI  block. 


Figure  13.  Hard-wired  data  cables 
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The  NI  block  is  moved  with  the  dual  canister  apparatus  each  time  it  is  removed 
from  the  carriage,  and  connects  to  the  chassis  (NI  PXIe-1082)  via  a  National 
Instruments  cable. 

2. 1.4.2.  Force  sensors  power  supply 

The  force  sensors  had  previously  been  powered  via  the  left  hand  side  of  a 
power  supply  on  the  carriage,  a  BK  precision  1673  triple  output  DC  power 
supply.  This  required  the  operator  to  twist  together  the  solid  orange  and 
green/white  wires  from  each  force  sensor  data  cable  to  clip  to  the  red  alligator 
clip,  leading  to  the  power  source.  The  black  alligator  clip  from  the  power  source 
would  then  clip  to  black  wires  that  needed  to  twist  together  from  the  NI  block.  If 
the  operator  wanted  to  re-zero  the  force  sensors,  the  left  hand  side  voltage  dial  of 
the  BK  power  supply  was  simply  turned  to  zero  to  power  down  the  force  sensors, 
and  then  turned  back  up  to  20V.  The  new  method  of  providing  power  to  the  force 
sensors  is  from  the  chassis,  through  a  NI  PXI-4110  programmable  DC  power 
supply.  Wires  go  from  the  chassis  to  clip  to  the  NI  block,  and  then  power  is 
conducted  through  the  force  sensor  data  cables  to  the  force  sensors.  Now,  to  re¬ 
zero  the  force  sensors  (which  was  done  at  the  start  of  each  data  collection  run), 
the  operator  simply  flips  a  switch  that  is  incorporated  into  the  force  sensor 
collection  software  (Lab VIEW  12.0f3). 

2. 1.4. 3.  LDM  power  supply 

The  LDM  used  to  be  powered  from  the  chassis,  the  20V  programmable 
supply  powers  the  force  sensors,  as  just  described.  The  LDM  is  now  powered  via 
the  left  hand  side  of  the  BK  power  supply. 
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2. 1.4.4.  Dual  canister  power  supply 

Power  is  supplied  to  the  roll  canister  via  an  Impulse  waterproof 
connector.  These  cables  were  previously  clipped  using  alligator  clips  that  were 
plug  in  leads  to  the  right  hand  side  of  the  BK  power  supply.  The  plug  in  leads 
were  cut  off,  and  directly  soldered  to  the  ends  of  the  red  and  black  cable  pair. 

The  negative  supply  line  was  grounded  to  earth/  carriage,  resulting  in  a  dramatic 
decrease  in  noise  recorded  at  the  DAQ  card. 

2.1.5.  Calibration  Improvements 
2. 1.5.1.  Previous  calibration  procedure 

The  previous  method  of  calibration  is  described  in  (Rauworth,  2014), 
chapter  3.2  and  chapter  4.  The  method  used  an  extra  foil  that  had  six  notches  cut 
into  the  leading  edge  at  known  distances  away  from  the  roll  axis.  A  set  of  three 
weights  was  hung  from  each  notch.  The  following  shows  a  typical  calibration  in 
progress: 


Figure  14.  Weights  hung  from  set  of  six  notches  (Rauworth,  2014) 
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This  procedure  was  repeated  for  the  positive  and  negative  direction  for  all  three 
coordinate  directions.  The  output  was  18  data  points  for  each  coordinate 
direction.  For  example,  if  the  weight  was  hung  in  the  positive  z-direction,  the 
output  was  for  one  quadrant  of  a  chart,  with  a  small  response  in  x,  a  small 
response  in  y,  and  a  large  response  in  z.  The  large  response  in  z  gave  three 
distinct  trend  lines  for  the  three  weights,  but  those  trend  lines  when  taken  as  a 
whole,  gave  a  different  straight  line.  The  slope  of  that  line,  fit  through  all  18 
points,  was  found  to  be  the  sensitivity  of  z  for  that  sensor  due  to  an  applied  z- 
force.  The  cross-terms  (sensitivity  of  x  and  y  due  to  an  applied  z-force)  were 
small  and  neglected.  In  summary,  this  procedure  used  a  force  applied  at  a 
moment  arm  to  produce  a  voltage  response  at  each  sensor.  Because  weights  and 
moment  arm  distances  were  all  known,  each  voltage  response  could  be  equated 
to  a  force  that  produced  it,  through  a  linear  relationship. 

2. 1.5. 2.  New  calibration  procedure 

Instead  of  using  three  weights  to  apply  various  moments  giving  sensor 
responses,  this  work  used  seven  weights  to  apply  force  in  each  coordinate  axis 
direction,  by  hanging  the  weights  directly  from  each  igubal  bearing  (since  the 
bearings  were  secured  above  each  force  sensor).  The  entire  calibration  section  of 
this  work  is  written  from  the  perspective  of  the  sensor  frames  of  reference, 
unless  specified,  which  are  shown  below.  The  sensor  directions  were  defined  as 
shown,  notice  that  they  do  not  adhere  to  a  right  hand  rule  coordinate  system. 
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Figure  15.  Sensor  coordinate  frames 

As  mentioned  earlier  it  was  prudent  to  increase  the  separation  between  them 
and  the  new  installation  method  required  the  sensors  to  be  installed  as  shown  in 
order  to  accommodate  the  pitch  motor.  For  reference,  the  body  frame  coordinate 
system  is  shown  in  blue.  Sensor  1,  on  the  port  side  of  the  dual  canister  and  closer 
to  the  body  frame  origin,  is  shown  in  red.  Sensor  2,  on  the  starboard  side  of  the 
dual  canister  and  closer  to  the  foil,  is  shown  in  green.  For  calibration,  the  x-  and 
z-coordinate  directions  were  clean  and  easy  to  obtain  calibration  data  from.  For 
x,  the  weights  were  hung  directly  out  the  side  hole  where  the  foil  would  be. 
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Figure  16.  Applying  an  x-direction  force 


For  z,  the  weights  were  hung  directly  out  the  top  of  the  delrin  housing. 
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Figure  17.  Applying  a  z-direction  force 


For  y,  it  was  more  complicated  because  the  delrin  housing  of  the  pitch  canister 
would  not  allow  any  direct  y-force  to  be  applied.  To  obtain  calibration  for  y,  the 
dual  canister  had  to  be  propped  up  at  a  45°  to  the  horizontal,  and  the  weights 
were  hung  obtaining  a  sensor  response  which  was  due  to  simultaneously  applied 
y-  and  z-forces.  The  procedure  to  find  the  responses  due  to  only  y-forces  is 
described  in  following  sections. 
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Figure  18.  Applying  simultaneous  y-  and  z-direction  forces 

Additionally,  the  new  calibration  procedure  does  not  neglect  the  cross-terms, 
however  small,  and  their  effect  is  used  through  all  calculations  to  find  lift  and 
thrust  forces  on  the  flapping  foil.  The  only  disadvantage  to  directly  hanging 
weights  from  the  bearings  is  that  only  one  direction  of  applied  force  could  be 
used  to  find  the  sensitivity  matrices  for  each  sensor.  For  sensor  axis  x,  force 
could  only  be  applied  out  the  foil  hole,  thus  it  was  a  negative  y-force  in  the  body 
coordinate  frame.  For  sensor  axis  z,  force  could  only  be  applied  out  the  open 
delrin  top  of  the  pitch  canister,  thus  it  was  a  negative  z-force  in  the  body 
coordinate  frame.  For  sensor  axis  y,  it  was  in  the  forward  direction  out  the  open 
delrin  top  of  the  pitch  canister,  thus  it  was  a  positive  x-force  in  the  body 
coordinate  frame. 

2. 1.5.3.  Sensitivity  and  calibration  matrices 
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For  both  the  sensitivity  and  calibration  matrices,  the  following 
nomenclature  is  used: 

Ab^ 

™ cd 

Equation  8.  Example  of  nomenclature 

"A"  can  be  either  "Se”  or  "Ca"  to  identify  itself  as  an  element  in  a  sensitivity  or  a 
calibration  matrix,  respectively,  "b"  can  either  be  "1"  or  "2”  to  identify  whether 
the  first  or  second  force  sensor  is  being  referred  to.  "c"  can  either  be  "x",  "y",  or 
"z"  to  identify  the  direction,  in  the  particular  sensor  coordinate  frame,  of  the 
sensor  response.  Lastly,  "d"  can  be  either  "x",  "y",  "z",  or  "yz"  to  identify  the 
direction,  in  the  particular  sensor  coordinate  frame,  of  the  applied  force.  Putting 
that  all  together,  the  following  element: 

Equation  9.  Nomenclature  in  use 

reads  as  sensor  1  sensitivity  in  the  z-direction  as  a  result  of  an  applied  x-force.  As 
stated  earlier,  the  sensitivity  matrix  elements  were  obtained  from  the 
appropriate  slopes.  Using  the  example  above  in  Equation  9,  that  particular 
element  will  be  obtained  from  the  figure  directly  below.  In  Figure  19,  the  left 
chart  for  sensor  1  has  a  red  linear  trend  line  which  corresponds  to  the  sensitivity 
element  in  question.  Units  of  sensitivity  are  Volts/Newtons. 
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Voltage,  V  Voltage,  V  Voltage, 


SI  response  to  applied  force  on  X 


S2  response  to  applied  force  on  X 


Figure  19.  Sensor  responses  due  to  applied  x-force 


SI  response  to  applied  force  on  Z 


S2  response  to  applied  force  on  Z 


Figure  20.  Sensor  responses  due  to  applied  z-force 


SI  response  to  applied  force  on  Y  and  Z  S2  response  to  applied  force  on  Y  and  Z 


Figure  21.  Sensor  responses  due  to  applied  y-  and  z-forces 


The  raw,  uncorrected  sensitivity  matrices  were  found  to  be: 
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Equation  10.  Raw  sensitivity  matrices 


The  middle  columns  are  sensitivities  due  to  the  simultaneously  applied  y-  and  z- 
force,  as  discussed  above.  For  sensor  1,  the  sensitivity  due  to  only  y-applied  force 
was  found  to  be  the  sensitivity  due  to  applied  yz-force  minus  the  sensitivity  due 
to  applied  z-force.  They  are  subtracted  because  in  the  sensor  1  coordinate  frame, 
negative  z-force  gives  a  negative  y-response. 
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Equation  11.  Sensor  1  correction  for  applied  y-force 


Similarly,  for  sensor  2  the  sensitivity  due  to  only  y-applied  force  was  found  to  be 
the  sensitivity  due  to  applied  yz-force  plus  the  sensitivity  due  to  applied  z-force. 
This  time,  in  the  sensor  2  coordinate  frame,  negative  z-force  gives  a  positive  y- 


response. 
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Equation  12.  Sensor  2  correction  for  applied  y-force 
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Thus  the  following  has  been  obtained  for  the  corrected  sensitivity  matrices. 
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Equation  13.  Final  sensitivity  matrices 


Allowing  the  following  response  plots  to  be  produced. 

SI  response  to  applied  force  on  Y  S2  response  to  applied  force  on  Y 


Figure  22.  Sensor  responses  due  to  applied  y-force 


What  remains  now  is  to  obtain  the  calibration  matrices  from  those  sensitivity 
matrices. 
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Equation  14.  Calibration  matrices 
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The  calibration  matrices  are  simply  the  inverse  of  the  sensitivity  matrices.  Units 


of  calibration  are  Newtons/Volts.  The  calibration  matrices  are  used  together 
with  the  output  voltages  of  each  sensor  in  order  to  find  the  output  forces 
experienced  at  each  sensor.  The  final  step  is  to  multiply  the  calibration  matrix  by 
the  output  voltages: 
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Equation  15.  Obtaining  forces  from  calibration  matrices  and  sensor  output  voltages 


Thus  each  force  output  in  a  given  direction  contains  contributions  from  all  three 
directions,  since  there  is  cross-sensitivity  seen  in  the  sensors.  The  equation 
above  is  listed  for  sensor  1  and  the  same  equation  is  used  for  sensor  2.  Note  the 
output  force  directions  are  in  the  individual  sensor  reference  frame,  the 
directions  were  shown  above  in  Figure  15.  This  means  that  there  will  still  be  a 
transformation  required  from  sensor  frame  forces  to  get  body  frame  forces.  The 
transformation  is  described  later. 

Although  the  method  of  calibration  was  improved  in  this  work,  errors  may 
still  have  been  introduced.  The  force  sensors  require  in  situ  calibration  once 
installed,  since  they  do  not  arrive  calibrated  from  the  manufacturer.  The  new 
calibration  method  directly  applies  force  to  each  bearing,  attached  above  each 
force  sensor.  To  accomplish  that  though,  the  pitch  shaft  had  to  be  removed.  To 
re-install  the  pitch  shaft,  the  bearings  needed  to  be  unscrewed  from  their 
mounting  plates  above  the  force  sensors,  slipped  over  the  pitch  shaft,  and  then 
placed  back  into  the  delrin  housing.  Errors  were  possibly  introduced  due  to  a 
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change  in  installation  after  calibration.  Additionally,  only  one  direction  of  each 
coordinate  axis  could  be  calibrated,  since  the  physical  arrangement  of  the  pitch 
canister  did  not  allow  force  to  be  applied  in  the  opposite  direction.  The 
assumption  was  that  the  sensor  response  was  linear  in  the  opposite  direction. 

2.2.  Experimental  Method 
2.2.1.  Experimental  Setup 

The  large  tow  tank  in  the  Sheets  Laboratory  on  the  Bay  Campus  was  used 
for  this  work,  the  most  important  aspect  of  which  was  to  configure  the  'beach'.  A 
wave  generator  (not  used  for  this  work)  is  installed  at  the  beginning  of  this  30m 
tow  tank,  and  the  floor  of  the  tow  tank  is  called  a  beach  because  it  is  often 
configured  as  a  sloping  beach,  to  observe  how  waves  run  up  or  break.  The  tow 
tank  is  actually  a  rectangular  prism,  cast  in  concrete,  and  the  beach  is  comprised 
of  metal  plates.  Each  metal  plate  is  2.5m  long  and  there  are  seven  plates  at  the 
end  of  the  tow  tank.  The  metal  plates  are  configured  by:  attaching  winches  to  the 
sides  of  the  tank,  clipping  the  winch  strap  end  into  the  seam  between  each  metal 
plate,  tightening  the  winch  strap  as  a  pair  to  ensure  uniformity,  unscrewing  the 
large  stainless  steel  bolts  that  pin  the  metal  plate  seams  into  the  inside  concrete 
wall  of  the  tow  tank,  using  the  winch  to  raise  or  lower  the  seam  as  appropriate, 
screwing  the  stainless  steel  bolts  back  in  to  pin  the  metal  plate,  and  then 
loosening  and  removing  the  winch  strap. 
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Figure  23.  Winches  ready  to  move  the  metal  plate  seams  of  the  beach 


The  tow  tank  was  configured  with  the  profile  of  water  depths  as  depicted  below: 


This  gave  a  deep  portion  of  8m  of  run  (tank  bottom  and  plate  1),  a  transition 
portion  of  2.5m  (plate  2),  and  then  a  shallow  portion  of  7m  of  run  (plates  3,  4,  and 
5).  Because  the  carriage  was  run  from  deep  to  shallow,  in  one  continuous  run  the 
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flapping  foil  was  able  to  experience  flapping  in  the  middle  of  the  water  column 
(hereafter  referred  to  as  "freestream”)  and  flapping  very  close  to  the  beach 
("near  bottom”).  For  the  freestream  case  at  the  roll  axis  the  H/c=8.0,  indeed,  the 
water  depth  was  so  deep  relative  to  the  foil  chord  that  if  the  foil  was  rotated  to 
point  straight  down  toward  the  beach,  the  tip  would  still  be  at  H/c=2.8.  Thus  for 
the  freestream  case,  there  is  absolutely  no  influence  from  the  ground,  nor  from 
the  walls  of  the  tow  tank.  In  the  freestream  case,  the  foil  flaps  free  and  clear  in 
deep  water.  When  the  foil  is  at  a  roll  amplitude  of  zero,  it  is  0.32m  below  the 
water  free  surface  (an  H/c  equivalent  of  3.2  below  water). 

In  order  to  get  the  largest  magnitude  response  for  ground  effect,  the 
desire  was  to  have  the  dual  canister  move  as  close  to  the  beach  as  possible.  Due 
to  the  radius  of  the  dual  canister  itself  and  the  metal  support  bracket,  it  is  only 
possible  to  get  the  roll  axis  in  ground  effect  to  a  minimum  H/c=l.l,  that  is  a 
physical  limitation.  It  is  a  realistic  limitation  however,  since  the  same  and  even 
larger  distances  away  from  the  ground  would  be  encountered  on  deployed  UUVs 
that  use  flapping  foils  for  propulsion,  due  to  the  associated  mechanisms 
necessary  for  roll  and  pitch.  Using  a  combination  of  the  fixed  length  carriage 
attachment  structure  and  wooden  blocks  to  shore  it  up  on  the  tow  tank  carriage 
itself,  the  metal  support  bracket  of  the  dual  canister  was  able  to  be  placed  1cm 
away  from  the  ground.  Note  that  the  overall  length  of  the  carriage  attachment 
with  dual  canister  attached  was  1.38m,  and  its  weight  was  in  excess  of 
approximately  350N.  Note  the  following  photograph,  showing  a  front  view  of  the 
dual  canister  at  the  closest  configuration  to  the  ground  (looking  in  the  body 
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frame  system  in  negative  x-direction).  The  fin  tip  is  touching  the  ground,  at  a  roll 


amplitude  of  a  little  over  13°. 


Figure  25.  Dual  canister  closest  to  the  beach 


Figure  26.  Side  view  of  dual  canister  near  bottom 


36 


The  maximum  roll  amplitude  can  only  be  12°,  and  that  value  has  been  used  for  all 
experiments  in  this  work.  Care  was  taken  to  ensure  that  the  bottom  of  the  dual 
canister  would  not  impact  the  beach.  The  most  important  aspect  was  the  safety 
and  security  of  school  assets  and  equipment.  To  establish  the  height  of  wooden 
blocks  needed  in  order  to  get  the  foil  as  close  as  possible  to  the  beach,  the 
carriage  was  operated  by  one  individual  outside  the  tow  tank.  The  author  was 
submerged  in  the  tow  tank,  next  to  the  dual  canister,  while  the  carriage  went 
along  at  a  snail's  pace.  Height  of  blocks  and  final  metal  plate  configuration  in  the 
tank  was  determined  such  that  the  dual  canister  will  never  impact  the  beach,  as 
long  as  the  carriage  is  stopped  above  metal  plate  5.  Foil  roll  amplitude  was 
limited  so  that  only  fluid  forces  over  the  foil  were  measured  and  not  forces  due  to 
the  fin  tip  impacting  the  beach.  With  roll  amplitude  fixed,  the  heave  to  chord 
ratio  is  also  therefore  fixed,  found  to  be  ho.7/c=0.91  for  all  experiments. 

2.2.2.  Experimental  Procedure 

The  procedure  for  data  collection  runs  was  identical  to  the  procedure 
presented  in  (Rauworth,  2014),  with  one  small  change.  In  this  work,  the  force 
sensors  were  powered  down  and  back  up  again  prior  to  each  data  run,  using  the 
switch  on  the  force  sensor  collection  software. 

2.2.3.  Data  Processing 

Many  of  the  same  data  processing  techniques  from  (Rauworth,  2014) 
were  applied  in  this  work.  As  described  above,  the  results  from  calibration 
produce  calibration  matrices  which  are  multiplied  by  sensor  voltages  to  give 
forces  experienced  at  each  sensor.  Only  once  sensor  forces  have  been  found  are 
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they  translated  into  the  body  frame  coordinate  system.  Noting  all  respective 


directions  from  Figure  15,  it  is  easy  to  see  how  the  following  translation  is  made 
from  sensor  forces  to  body  forces,  for  each  respective  sensor: 
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Equation  16.  Translation  from  sensor  frame  to  body  frame 


The  flapping  foil  generates  lift  and  thrust  forces  as  it  flaps.  The  lift  force  is  taken 
to  act  through  the  assumed  hydrodynamic  center,  some  distance  away  from  the 
origin.  That  distance  is  actually  ro.7,  given  in  Equation  2.  Lift  applies  a  moment 
about  the  body  x-axis.  That  moment  is  known,  because  the  force  sensors  have 
recorded  the  body  forces,  and  each  sensor  is  a  known  distance  away  from  the 
origin.  Thus  the  trail  to  find  lift  force  at  the  foil  is  as  follows:  record  forces 
encountered  at  sensor,  find  the  moment  at  the  origin  that  would  produce  those 
forces,  divide  out  ro.7,  and  finally  obtain  lift  force  at  the  hydrodynamic  center. 

The  same  procedure  is  used  to  find  the  thrust  force  developed  at  the  foil  assumed 
hydrodynamic  center,  which  will  produce  a  moment  about  the  body  z-axis.  The 
following  representation  has  not  been  drawn  to  scale.  It  shows  forces  and 
dimensions  related  to  the  pitch  canister  only  illustrating  how  moments  can  be 
calculated  to  the  origin,  and  then  the  lift  and  thrust  forces  can  be  backed  out.  The 
number  superscripts  refer  to  either  sensor  1,  or  sensor  2,  depending  on  which 
corresponding  number  is  used. 
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front 


Figure  27.  Three  view  drawing  to  represent  pitch  canister 


From  that  figure,  it  is  easy  to  write  out  the  calculations  which  must  be  performed 

to  find  lift  and  thrust  forces: 

r  _  ^(-z1)  +  Fy2(-z2)  +  F/Cy1)  +  Fz2(-y2) 

" 


T  ^(-y1)  +^2(y2) 

r0.7 

Equation  17.  Equations  to  find  lift  and  thrust  forces 

At  this  point  it  is  important  to  mention  that  the  largest  assumption,  and  hence 
potential  source  of  error  in  this  work,  was  in  the  assumed  location  of  the 
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hydrodynamic  center  of  the  foil.  The  same  location,  ro.7  has  been  consistently 
used  in  all  the  preceding  works  that  have  been  built  upon  here,  such  as  in 
(Polidoro,  2003),  (Techet,  2008),  and  (Rauworth,  2014).  The  location  of  the 
actual  hydrodynamic  center  was  still  not  able  to  be  accurately  determined,  due  to 
the  low  resolution  and  high  noise  still  present  in  the  system.  Looking  at  the 
equations  above  it  is  easy  to  see  that  lift  and  thrust  forces  calculated  can  change 
greatly  due  to  the  denominator.  The  maximum  forces  that  might  be  calculated 
would  be  found  if  the  denominator  used  was  the  radius  at  the  root  of  the  foil, 
similarly,  the  forces  would  be  minimum  if  the  radius  at  the  tip  of  the  foil  was 
used.  That  fact  provides  two  bounds,  that  the  lift  and  thrust  forces  are  to  be 
found  in  between. 

As  mentioned  before,  one  run  down  the  tow  tank  consists  of  three  distinct 
flapping  regimes;  namely  freestream,  transition  zone,  and  near  bottom.  It  was 
very  important  to  synchronize  the  LDM  readout  with  the  data  collected  from  the 
force  sensors  for  the  next  step  in  data  processing.  Once  the  lift  and  thrust  forces 
are  known  for  the  entire  run  they  were  phase  averaged  in  order  to  produce  plots 
that  showed  the  average  forces  over  one  flapping  cycle.  Phase  averaging  was 
conducted  twice  for  each  run,  first  for  flapping  in  the  freestream,  and  second  for 
flapping  near  bottom.  The  processing  code  also  calculated  the  maximum  values 
and  minimum  values  of  lift  and  thrust  force,  as  well  as  the  mean  lift  and  mean 
thrust  coefficients  for  each  run,  for  both  freestream  and  near  bottom.  All  of  these 
values  will  be  used  in  the  next  chapter  to  quantify  the  ground  effect  that  was 
observed. 
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3.  Results  and  Discussion 


3.1.  Establishing  zero  mean  lift  coefficient 

Prior  to  collecting  data  to  process,  a  set  of  runs  were  conducted  in  order  to 
determine  the  pitch  bias  that  would  give  a  zero  mean  lift  coefficient.  The  zeroing 
method  for  this  flapping  foil  has  always  proven  problematic  because  the  setup 
inside  the  pitch  canister  is  not  perfect  and  allows  for  a  few  degrees  of  pitch 
change.  The  belt  that  drives  the  foil  to  pitch  has  been  installed  as  tightly  as 
possible,  nonetheless,  when  the  dual  canister  is  powered  up  and  supposedly 
locked  in  position  it  is  still  possible  to  rotate  the  foil  in  pitch  by  approximately  a 
degree.  A  set  of  11  runs  were  conducted  for  0bias=-lOo:2°:lO°  at  St=0.5  and 
ocmax=30°.  The  goal  was  to  determine  the  pitch  bias  for  which  the  flapping  foil 
would  exhibit  zero  mean  lift  coefficient  during  the  near  bottom  portion  of  the 
run.  To  the  nearest  integer  it  was  found  to  be  3°  pitch  bias  and  was  set  as  the 
new  zero  pitch  bias: 
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Figure  28.  Mean  lift  coefficient  as  a  function  of  pitch  bias 
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Note  that  to  find  the  pitch  bias  at  which  the  mean  lift  coefficient  was  zero,  only 
one  combination  of  St  and  amax  was  used.  This  means  that  the  pitch  bias  was 
accurate  only  for  that  kinematic  condition.  This  may  have  introduced  some  error 
in  that  a  3°  pitch  bias  was  used  as  zero,  to  represent  all  data  points  when  flapping 
near  bottom.  It  would  have  been  prohibitive  in  terms  of  time,  though,  to  run 
through  a  pitch  bias  sweep  for  each  of  the  combinations  of  St  and  amax  to  find  a 
zero  mean  lift  coefficient  pitch  bias  for  each.  Additionally,  the  nearest  integer 
pitch  bias  was  used,  but  it  should  have  been  to  a  fraction  of  a  degree.  A  positive 
pitch  bias  is  defined  as  the  pitch  amplitude  centerline  being  above  the  horizontal 
at  the  leading  edge,  so  that  at  zero  pitch  amplitude  the  foil  has  some  positive 
angle  of  attack.  This  also  means  that  the  pitch  centerline  at  the  trailing  edge  is 
pointed  toward  the  ground  at  a  positive  pitch  bias.  Figure  28  above  shows  the 
variation  in  mean  lift  coefficient  with  close  approach  to  ground.  At  zero  pitch 
bias,  flapping  near  bottom  has  a  slightly  positive  lift  coefficient,  Cl=0.022, 
whereas  flapping  in  the  freestream  has  a  negative  lift  coefficient,  Cl=-0.087.  The 
result  in  Figure  28  agrees  with  the  observation  in  (Licht  &  Dahl,  2013)  that  the 
freestream  case  must  have  a  larger  than  1°  positive  pitch  bias  in  order  to 
generate  the  same  lift  as  the  near  bottom  case.  This  effect  magnifies  as  larger 
mean  lift  coefficients  are  sought,  such  that  at  a  mean  Cl=0.4,  the  difference  is  a  2° 
positive  pitch  bias. 

3.2.  Contour  plots 

The  experimental  matrix  of  tests  was  conducted  immediately  after  finding 
the  pitch  bias  zero,  without  changing  the  roll  or  pitch  position  of  the  dual 
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canister.  Each  run  was  conducted  at  the  pitch  bias  value  mentioned  above,  which 


would  nominally  give  zero  mean  lift  coefficient  for  the  near  bottom  portion.  For 
reference  the  following  is  a  table  listing  the  flapping  frequency,  roll  amplitude 
(AMX),  pitch  amplitude  (AMY),  and  corresponding  Strouhal  number  and 
maximum  angle  of  attack.  The  row  that  is  bold  is  the  test  condition  at  which  the 
zero  mean  lift  coefficient  test  above  was  conducted  at. 


Test  No 

f  (Hz) 

AMX  (°) 

AMY  (°) 

St 

OCmax  C  ) 

1 

0.82 

12 

23 

0.3 

20 

2 

0.82 

12 

18 

0.3 

25 

3 

0.82 

12 

13 

0.3 

30 

4 

0.82 

12 

8 

0.3 

35 

5 

1.1 

12 

32 

0.4 

20 

6 

1.1 

12 

27 

0.4 

25 

7 

1.1 

12 

22 

0.4 

30 

8 

1.1 

12 

17 

0.4 

35 

9 

1.37 

12 

39 

0.5 

20 

10 

1.37 

12 

33 

0.5 

25 

11 

1.37 

12 

28 

0.5 

30 

12 

1.37 

12 

22 

0.5 

35 

13 

1.65 

12 

47 

0.6 

20 

14 

1.65 

12 

40 

0.6 

25 

15 

1.65 

12 

33 

0.6 

30 

16 

1.65 

12 

27 

0.6 

35 

Table  1.  Experimental  matrix 


Having  conducted  all  of  the  runs  listed  above,  the  following  contour  plots  were 
produced. 
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Figure  29.  Contour  plot  for  mean  lift  coefficient  in  freestream  (H/c=8.0) 


Figure  30.  Contour  plot  for  mean  lift  coefficient  near  bottom  (H/c=l.l) 
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Figure  31.  Contour  plot  for  change  in  mean  lift  coefficient 

Notice  immediately  that  in  Figure  29  almost  every  data  point  is  a  negative  mean 
lift  coefficient,  again  highlighting  that  in  the  freestream,  there  will  be  a  negative 
lift  force  generated.  The  only  region  for  which  flapping  in  the  freestream 
generates  positive  lift  is  for  the  highest  St  and  lowest  amax,  which  indicates  most 
aggressive  flapping  motion,  having  the  largest  frequency  and  the  largest  pitch 
amplitude.  In  Figure  30,  the  opposite  is  true,  meaning  that  there  is  only  a  small 
region  of  three  data  points  where  the  mean  lift  coefficient  is  negative.  Perhaps 
most  helpful  is  to  look  at  Figure  3lFigure  31,  which  shows  the  difference  between 
flapping  near  bottom  and  in  the  freestream.  The  data  is  plotted  for  Cl  near 
bottom  minus  Cl  in  freestream.  It  is  clear  that  at  every  point  in  the  test  matrix 
the  mean  lift  coefficient  near  the  bottom  is  higher  than  the  mean  lift  coefficient  in 
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the  freestream.  This  indicates  that  the  foil  is  in  ground  effect  for  every  case 
tested.  The  largest  change  in  mean  lift  coefficient  is  for  St=0.3  and  amax=20°,  and 
then  for  a  region  where  St=0.5  and  amax=30°  and  35°.  The  smallest  change  in 
mean  lift  coefficient  is  to  be  found  for  St=0.6  and  amax=20°.  Notice  that  for  a  given 
ocmax=20°,  the  largest  and  smallest  mean  lift  coefficient  changes  are  found.  For  the 
case  of  low  St  (low  frequency),  each  flap  is  long  in  duration,  producing  a  short  but 
broad  lift  trace  as  a  function  of  time.  For  the  case  of  high  St  (high  frequency), 
each  flap  is  short  in  duration,  producing  a  taller  but  narrower  lift  trace  as  a 
function  of  time.  Because  low  St  equates  to  larger  change  in  mean  lift  coefficient 
than  high  St,  the  value  for  mean  lift  must  be  larger  for  low  St,  due  to  that  longer 
duration  flap. 

The  boost  in  lift  due  to  ground  effect  has  been  revealed,  so  now  to 
determine  if  there  is  a  benefit  to  the  thrust  produced.  Based  on  aerospace 
engineering,  the  expectation  is  there  will  be  some  benefit,  as  induced  drag  is 
decreased  in  ground  effect.  However,  the  study  by  (Blevins  &  Lauder,  2013) 
showed  that  there  is  not  necessarily  a  locomotor  advantage.  Thrust  force 
produced  by  the  foil  on  the  dual  canister  has  been  defined  as  positive  forward, 
namely  positive  in  the  body  x-direction. 
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Figure  32.  Contour  plot  for  mean  thrust  coefficient  in  freestream  (H/c=8.0) 


Figure  33.  Contour  plot  for  mean  thrust  coefficient  near  bottom  (H/c=l.l) 
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Figure  34.  Contour  plot  for  change  in  mean  thrust  coefficient 

Looking  at  Figure  32  first,  the  mean  thrust  coefficient  is  positive  for  the  entire 
plot,  which  means  that  there  is  positive  thrust  being  produced.  The  same  is 
evident  in  Figure  33,  the  case  for  near  the  bottom.  Again,  what  is  most  beneficial 
is  to  see  the  change  in  mean  thrust  coefficient  from  near  bottom  to  freestream 
conditions,  Figure  34.  What  is  evident  from  the  plot  is  that  the  largest  thrust 
benefit  to  being  in  ground  effect  is  found  for  a  St=0.6  and  amax=30°.  Surrounding 
that  data  point  is  the  only  region  that  is  positive,  so  there  is  only  going  to  be  a 
thrust  benefit  at  high  St  and  high  amax,  which  equates  to  a  rapid  flapping 
frequency  and  mid  to  low  range  pitch  amplitude.  Notice  that  there  are  large 
regions  of  negative  change  in  mean  thrust  coefficient,  which  means  that  less 
thrust  is  produced  near  the  bottom  as  compared  to  flapping  in  the  freestream. 
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These  areas  are  for  all  conditions  where  St=0.3,  then  for  all  conditions  where 
ocmax=20°  and  25°,  then  again  for  a  St=0.4  and  amax=35°.  This  actually  means  that 
there  is  more  drag  force  near  bottom  on  the  foil,  than  there  was  in  the 
freestream.  Just  as  observed  in  (Blevins  &  Lauder,  2013),  they  found  that 
undulating  fins  generally  incur  costs  from  moving  close  to  a  solid  boundary. 


3.3.  Results  for  different  series  of  Strouhal  number 


Figure  35.  Change  in  mean  lift  coefficient  as  a  function  of  maximum  angle  of  attack 


The  change  in  mean  lift  coefficient  for  every  value  of  St  and  amax  in  the 
experimental  range  is  a  positive  number;  mean  lift  coefficient  near  bottom  is 
greater  than  mean  lift  coefficient  in  the  freestream.  For  an  amax  of  25°  to  35°,  the 
change  in  mean  lift  coefficient  is  in  a  fairly  narrow  band.  It  is  only  for  an 
ocmax=20°,  so  largest  pitch  amplitude  and  most  feathered  foil,  that  the  results 
diverge  the  most.  The  results  for  St=0.3  and  St=0.6  have  opposite  trends,  concave 
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up,  as  opposed  to  concave  down.  The  middle  two  St  have  fairly  stable  trends, 
such  that  there  does  not  appear  to  be  much  difference  to  changing  the  amax  over 
the  range  studied.  The  data  in  Figure  35  shows  what  was  observed  earlier  in 
Figure  31,  that  if  the  largest  change  in  mean  lift  coefficient  is  sought,  it  is  to  be 
found  for  St=0.3  and  amax  of  20°.  Similarly,  the  smallest  change  is  for  a  St=0.6  and 

ffmax  of  20°. 


Figure  36.  Change  in  maximum  lift  force  as  a  function  of  maximum  angle  of  attack 

The  figure  directly  above  gives  another  indication  that  for  all  experiments  the 
highest  single  lift  force  in  ground  effect  was  larger  than  the  highest  single  lift 
force  in  the  freestream,  proving  the  lift  benefit  of  ground  effect.  For  St=0.3  and 
amax=20°,  the  combination  of  both  the  largest  change  in  mean  lift  coefficient,  and 
the  largest  change  in  maximum  lift  force  is  observed.  The  percent  difference 
between  the  maximum  lift  force  near  bottom  and  in  the  freestream  at  that  point 
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is  45%.  This  seems  like  a  huge  advantage,  however,  from  a  practical  point  of 
view  the  value  of  this  observation  is  limited  since  the  foil  is  producing  net  drag  at 
this  kinematic  condition. 


Figure  37.  Change  in  minimum  lift  force  as  a  function  of  maximum  angle  of  attack 


The  magnitude  of  the  change  is  smaller  when  considering  the  minimum  observed 
lift  force,  but  again,  all  changes  observed  are  positive.  In  the  case  of  minimum  lift 
force,  the  positive  change  means  that  the  minimum  lift  force  due  flapping  near 
the  bottom  is  less  negative  than  that  due  to  the  freestream,  so  indicates  more 
beneficial  boost  to  lift. 
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Figure  38.  Change  in  mean  thrust  coefficient  as  a  function  of  maximum  angle  of  attack 


The  figure  above  clearly  illustrates  which  combinations  of  St  and  amax  will 
generate  a  thrust  benefit  due  to  flapping  in  ground  effect,  namely  the  region 
where  the  change  in  thrust  coefficient  is  positive.  The  same  is  shown  in  the 
contour  plot,  Figure  34.  If  the  desired  operational  conditions  of  a  UUV  is  to 
generate  more  thrust  due  to  flapping  near  the  ground,  as  compared  to  flapping  in 
the  freestream,  that  is  the  area  to  choose.  For  St=0.3,  more  thrust  is  never 
produced,  all  that  is  obtained  is  a  larger  amount  of  drag  for  this  low  frequency 
flapping.  For  St=0.4,  the  mid-range  amax  values  give  a  thrust  benefit.  For  the 
higher  St,  it  is  only  at  higher  amax  values  that  a  thrust  benefit  is  seen,  but  not  for 
the  highest  St=0.6. 
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Figure  39.  Change  in  maximum  thrust  force  as  a  function  of  maximum  angle  of  attack 


The  figure  above  is  perhaps  a  little  misleading  at  first  glance.  It  shows  in  all  cases 
that  the  maximum  thrust  force  for  flapping  near  the  bottom  is  less  than  the  value 
for  flapping  in  the  freestream.  This  situation  can  still  result  in  an  overall  positive 
mean  thrust  coefficient,  as  it  is  calculated  by  taking  the  mean  thrust  force  value 
over  the  entire  flapping  cycle,  from  Equation  7.  This  will  be  explored  later  by 
observing  the  phase  averaged  thrust  force  curves. 
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Figure  40.  Change  in  minimum  thrust  force  as  a  function  of  maximum  angle  of  attack 

The  figure  above  is  also  slightly  misleading,  in  the  same  way  as  Figure  39  was. 
The  interpretation  will  be  revisited  after  showing  phase  averaged  thrust  force 
curves. 

3.4.  Results  for  a  single  maximum  angle  of  attack 

As  stated  earlier,  the  procedure  to  find  zero  mean  lift  coefficient  by 
running  through  a  sequence  of  pitch  bias  values  was  conducted  at  St=0.5  and 
ocmax=30°.  In  the  data  set,  there  were  also  three  other  points  at  that  same  amax. 
These  four  data  points  will  be  used  now  for  a  comparison,  as  St  is  varied.  The 
results  can  be  followed  by  referring  to  the  contour  plots,  and  tracing  a  horizontal 
line  through  ffmax—  30°. 
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Figure  41.  Change  in  mean  lift  coefficient  or  mean  thrust  coefficient  as  a  function  of 
Strouhal  number  for  maximum  angle  of  attack  of  30° 


The  figure  above  shows  that  in  all  four  cases,  the  change  in  mean  lift  coefficient  is 
positive  and  has  a  maximum  for  St=0.5.  This  means  that  ground  effect  is 
producing  a  lift  benefit.  The  mean  thrust  coefficient  is  positive  for  the  three 
higher  St,  and  slightly  negative  for  St=0.3.  These  results  agree  with  the  previous 
numerical  study  for  insect  wing  flapping  in  air  (Wu,  Shu,  Zhao,  &  Yan,  2014).  For 
all  cases  the  mean  lift  coefficient  increased  in  ground  effect,  and  the  thrust 
coefficient  shows  the  same  pattern,  namely  that  thrust  coefficient  decreased  for 
low  St,  and  increased  for  high  St. 
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Figure  42.  Change  in  minimum  lift  force  or  maximum  lift  force  as  a  function  of  Strouhal 
number  for  maximum  angle  of  attack  of  30° 


The  figure  above  shows  the  lift  forces  for  the  four  data  points.  Notice  that  all 
points  are  positive,  showing  the  improvement  in  lift.  The  largest  changes  in  max 
lift  force  are  for  St=0.4  and  0.5. 


Figure  43.  Change  in  minimum  thrust  force  or  maximum  thrust  force  as  a  function  of 
Strouhal  number  for  maximum  angle  of  attack  of  30° 
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The  figure  above  shows  the  changes  in  thrust  forces,  that  again  may  be 
misleading,  but  these  are  simply  the  maximum  and  minimum  values.  The  next 
section  will  clear  up  and  elaborate. 

3.5.  Phase  averaged  lift  and  thrust  force  plots 

The  following  plots  have  been  produced  by  taking  many  flapping  cycles 
and  phase  averaging  them  into  one  upstoke  and  downstroke,  for  the  freestream 
and  near  bottom  regimes.  Four  examples  are  given  for  the  four  data  points 
where  amax=30°.  Dashed  lines  are  for  freestream,  solid  are  for  near  bottom. 


Figure  44.  Phase  averaged  lift  and  thrust  force  for  St=0.3  and  amax=30° 
(dashed=freestream,  solid=near  bottom) 
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The  figure  above  shows  that  the  lift  force  is  consistently  greater  in  the  case  of 
near  bottom  than  it  is  for  the  freestream.  Notice  that  the  opposite  is  true  for  the 
thrust  force,  the  freestream  has  a  larger  positive  maximum  value.  For  this  case, 
as  seen  in  Figure  41,  the  change  in  mean  thrust  coefficient  is  slightly  negative,  so 
there  is  more  thrust  produced  during  flapping  in  the  freestream,  as  compared  to 
near  bottom. 


Figure  45.  Phase  averaged  lift  and  thrust  force  for  St=0.4  and  amax=30° 
(dashed=freestream,  solid=near  bottom) 

Next  considering  the  figure  above,  again  the  maximum  lift  force  for  near  bottom 
is  larger  than  in  freestream,  but  the  maximum  thrust  force  is  smaller.  For  this 
case  of  St=0.4,  the  change  in  mean  thrust  coefficient  is  actually  positive,  and  will 
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be  for  the  next  two  St  that  follow.  This  perfectly  explains  the  misleading  data 


from  earlier  (Figure  39,  Figure  40,  Figure  43),  because  here  the  first  peak  in 
thrust  has  a  lot  more  area  underneath  it  for  the  near  bottom  case,  and  the  second 
peak  does  not  have  such  a  large  difference  in  area  underneath.  The  first  peak 
increase  more  than  makes  up  for  the  second  peak  decrease,  and  the  net  result  is  a 
larger  mean  thrust  coefficient,  as  calculated  from  Equation  7. 


Figure  46.  Phase  averaged  lift  and  thrust  force  for  St=0.5  and  amax=30° 
(dashed=freestream,  solid=near  bottom) 


For  the  figure  above,  all  of  the  same  attributes  from  the  previous  two  cases  are 
seen.  Notice  that  the  peaks  in  thrust  are  moving  closer  to  the  same  height  as  St  is 
increased.  The  figure  also  corresponds  to  the  same  kinematic  conditions  as  the 
lift  and  thrust  plot  taken  from  (Rauworth,  2014),  which  is  Figure  8  in  this  work. 
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Notice  the  trends  in  lift  and  thrust  are  very  similar  to  Figure  8,  though  here  the 
magnitude  of  the  traces  are  a  few  Newtons  of  force  lower.  It  is  difficult  to 
determine  which,  if  any,  of  the  two  examples  compared  is  the  correct,  true 
representation  of  forces  generated  by  the  foil,  especially  since  the  moment 
equation  procedure  led  to  those  forces.  The  following  is  simply  a  side  by  side 
comparison  showing  Figure  8  and  Figure  46  again: 


Figure  47.  Comparison  between  lift  and  thrust  plot  for  (Rauworth,  2014)  and  this  work 

What  is  important  is  just  to  note  that  the  trends  are  the  same,  and  that  in  the 
current  work,  an  influence  in  the  data  from  the  foil  being  near  bottom  was  seen. 
Phase  averaging  was  conducted  slightly  differently  for  both  plots,  so  the 
expectation  is  that  the  lift  and  thrust  traces  would  not  line  up  for  perfect  side  by 
side  comparison. 
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Figure  48.  Phase  averaged  lift  and  thrust  force  for  St=0.6  and  amax=30° 
(dashed=freestream,  solid=near  bottom) 


For  the  figure  above,  notice  that  the  near  bottom  lift  trace  is  symmetrical.  For 
this  particular  case,  it  happens  to  be  that  the  largest  change  in  mean  thrust 
coefficient  is  observed  in  the  entire  data  set.  The  percent  difference  between 
maximum  thrust  force  near  bottom  and  in  the  freestream  is  6%,  so  if  an  overall 
thrust  benefit  due  to  ground  effect  is  sought,  that  change  in  instantaneous 
maximum  thrust  force  will  need  to  be  sensed.  The  percent  difference  in  mean 
thrust  coefficient  here  was  8%.  Of  course,  this  may  not  be  as  useful  as  it  seems, 
because  the  propulsive  efficiency  difference  has  not  been  calculated.  Similar  to 
findings  in  (Blevins  &  Lauder,  2013),  the  power  requirement  may  increase  at  this 
kinematic  condition  when  flapping  near  the  bottom. 
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4.  Design  Improvements  and  Error  Analysis 

4.1.  Design  Improvements 

The  experiments  indicate  that  despite  improvements  in  signal  to  noise 
ratio  over  previous  apparatus,  it  will  still  be  necessary  to  change  the  design 
inside  the  pitch  canister  to  achieve  quantifiable,  repeatable  results. 

4.1.1.  One  force  sensor 

The  most  ideal  method  of  force  sensing  would  be  to  only  have  one  sensor, 
instead  of  the  two  that  are  currently  installed.  One  larger  bearing  and  sensor 
could  be  installed  to  support  the  pitch  shaft,  and  there  would  no  longer  be 
confusion  about  which  direction  forces  were  acting  in,  and  having  to  resolve  both 
force  sensor  forces  into  one.  Also,  in  the  x-  and  y-directions  the  sensors  had  the 
same  range,  but  a  different  range  of  sensing  existed  for  the  z-direction.  Forces  in 
the  z-direction  were  very  important,  as  they  lead  to  the  lift  calculations  for  the 
foil. 

4.1.2.  Install  an  accurate  torque  sensor 

A  torque  sensor  should  be  installed  on  the  pitch  shaft.  With  the  present 
setup,  there  is  a  measure  of  torque  output  based  on  motor  current  for  the  roll  and 
pitch  motors,  but  the  data  rate  and  resolution  is  too  low  to  be  useful.  A  new 
torque  sensor  will  give  a  measure  of  the  power  required  to  drive  the  foil.  From 
that  data  the  propulsive  efficiency  can  be  determined  to  see  if  flapping  near  the 
bottom  is  more  efficient.  Those  kinematic  conditions  for  which  a  propulsive 
efficiency  boost  due  to  flapping  near  the  bottom  can  be  investigated  and 
employed  on  future  UUVs. 
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4.2.  Sources  of  error 


4.2.1.  Method  for  centering  foil  in  roll 

The  foil  was  pulled  up  to  the  support  bracket  limit,  and  then  commanded 
to  rotate  5200  counts  down  in  roll,  in  order  to  start  it  at  the  horizontal.  This 
method,  used  in  (Rauworth,  2014),  likely  introduced  some  error. 

4.2.2.  Pitch  bias  sweep  to  find  zero  mean  lift  coefficient 

In  future  tests,  a  zero  mean  pitch  bias  should  be  found  for  flapping  in  the 
freestream,  rather  than  near  bottom.  That  value  for  pitch  bias  should  also  be 
used,  to  two  decimal  places  rather  than  closest  integer,  for  all  experiments.  By 
doing  that,  the  "normal"  condition  would  be  flapping  in  the  freestream,  and  then 
any  ground  effect  by  flapping  near  the  bottom  would  be  seen  clearly  in  the  data. 

4.2.3.  Carriage  drive  mechanism  and  dual  canister  motors 

Another  source  of  error  may  be  found  in  the  carriage  drive  mechanism, 
which  is  old  and  fairly  unreliable.  The  LDM  should  be  connected  to  the  drive 
mechanism,  in  order  to  command  the  carriage  to  move  at  precise  speeds,  rather 
than  just  using  a  rheostat. 

4.2.4.  Motor  control 

Foil  position  has  been  plotted  based  on  the  input  command,  so  the  desired 
position  was  used.  As  mentioned  earlier,  the  roll  and  pitch  motors  of  the  dual 
canister  are  old,  and  their  actual  output  position  data  was  not  used.  A  potential 
source  of  error  is  therefore  introduced  if  the  desired  and  actual  motor  positions 
are  different. 
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5.  Future  Work 


5.1.  Stationary  carriage 

The  testing  platform  and  processing  software  was  only  setup  to  run  at 
some  non-zero  forward  speed.  It  would  be  valuable  to  determine  the  flapping 
conditions  required  for  station  keeping,  such  that  a  high  lift  is  required,  and 
basically  zero  thrust.  Station  keeping  for  a  UUV  may  be  employed  during 
swimmer  defense,  explosive  ordinance  disposal,  video  surveys  of  one  specific 
point  of  interest,  and  many  other  missions. 

5.2.  Under  a  free  water  surface 

Just  as  a  lift  and  thrust  benefit  was  observed  close  to  a  hard,  bottom 
surface,  an  investigation  could  be  conducted  into  the  effects  of  flapping  near  a 
free  surface.  Although  not  hard,  there  is  a  drastic  density  change  between  water 
and  air,  so  it  would  be  a  worthy  pursuit.  However,  the  wave  making  that  is  the 
result  of  flapping  near  the  free  surface  might  make  potential  gains  obsolete,  by 
wasting  energy.  The  damping  force  that  would  be  produced  would  dissipate 
energy. 

5.3.  Under  wave  conditions 

Real  world  deployment  of  UUVs  is  of  ultimate  importance  for  all  this 
investigation  work,  and  in  the  real  world  there  are  waves  to  deal  with  in  any 
significant  body  of  water.  The  tow  tank  that  testing  was  conducted  in  has  a  wave 
generator,  and  the  foil  performance  could  be  tested  under  diverse  wave 
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conditions  to  see  if  there  are  particular  wave  conditions  which  provide  benefits 
to  various  aspects  of  a  UUV  mission. 

5.4.  Near  walls 

Similar  to  flapping  near  a  horizontal  bottom  ground,  flapping  near  vertical 
walls  or  even  under  hard  surfaces  should  be  explored.  The  flapping  could  be 
oriented  to  have  a  flap  centerline  to  match  the  profile  of  the  hard  surface,  or  to 
flap  with  the  tip  of  the  foil  perpendicular  to  the  hard  surface  in  question.  The  two 
different  configurations  would  allow  either  most  of  the  span  to  interact  with  the 
hard  surface,  or  just  the  tip  so  that  end  effects  could  be  investigated.  Flapping 
perpendicular  to  a  hard  surface  may  provide  similar  benefits  to  comparing  three 
dimensional  flow  to  two  dimensional  flow  over  a  wing.  Tip  losses  are  minimized 
by  vertical  winglets  at  the  wing  tips  in  aerospace  applications,  a  similar  benefit 
remains  to  be  investigated  underwater. 

5.5.  Performance  and  operational  considerations 

In  future  work  many  more  parameters  to  find  the  kinematic  conditions  in 
which  a  benefit  from  ground  effect  is  observed  must  be  used,  such  as  power 
consumption  and  propulsive  efficiency.  A  much  more  comprehensive 
experimental  matrix  can  be  employed  to  get  more  resolution  by  varying  St  and 
amax  in  finer  increments.  Again  using  nature  as  an  inspiration,  observations  from 
flying  animals  utilizing  ground  effect  can  lead  to  some  operating  rules  for  UUVs  to 
maximize  the  benefit.  Some  animal  behaviors  have  been  identified  as:  (i)  an 
animal  cannot  fly  very  slowly  in  ground  effect;  (ii)  for  power  or  range  economy 
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an  animal  should  reduce  speed  when  commuting  in  ground  effect;  (iii)  an  animal 
should  show  less  tendency  to  vary  flight  speed  at  any  given  height  when  in 
ground  effect,  particularly  when  it  is  flying  slowly,  at  or  around  the  minimum 
power  speed;  and  (iv)  when  varying  flight  speed  over  a  surface  an  animal  should 
reduce  height  at  higher  flight  speeds  (Rayner,  1991).  The  same  paper  discussed 
even  more  observations  from  nature.  Pelicans  gliding  over  water  fly  at  a  much 
lower  altitude  (H/s=0.31)  than  when  they  are  flapping  (H/s=0.50).  Notice  that 
the  author  used  span  to  non  dimensionalize  the  altitude.  Skimmers  and  myotid 
bats  adopt  an  asymmetric  wingbeat  in  which  the  downwards  excursion  of  the 
wing  below  the  body  has  a  much  lower  amplitude  than  the  upper  half  of  the 
wingbeat.  In  this  work  an  equal  upstroke  and  downstroke  roll  amplitude  was 
used,  but  future  work  can  certainly  investigate  increasing  the  roll  amplitude  for 
the  foil  during  the  upper  half  of  the  flap,  when  it  is  away  from  the  ground. 

As  far  as  animals  swimming  in  water  near  the  bottom,  some  observations 
have  been  made  by  (Blevins  &  Lauder,  2013).  They  noted  that  live  animals 
change  their  behavior  when  swimming  close  to  a  solid  boundary.  Due  to  the 
stiffness  of  the  fin  model  used,  the  authors  were  not  able  to  detect  subtle 
influences  that  ground  effect  has  on  fin  shape.  They  noted  that  kinematic 
changes  observed  in  fish  swimming  in  ground  effect  are  due  to  active  modulation, 
not  passive  effects.  Furthermore,  these  changes  consisted  of  reduced  tailbeat 
frequency  and  amplitude,  which  may  be  a  mechanism  for  transforming  potential 
locomotor  penalties  into  performance  benefits.  Stingrays  have  extremely  fine 
control  of  pectoral  fin  conformation  with  the  potential  for  precise  kinematic 
tuning  to  exploit  ground  effect. 
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These  behavior  changes  noted  in  animals  exploiting  ground  effect  may 
provide  a  starting  point  for  more  investigation  into  this  topic.  Of  course  these 
behaviors  will  have  to  be  programmed  into  the  UUV  profile  in  order  to  take 
maximum  advantage  of  ground  effect.  When  the  UUV  is  required  to  go  into 
ground  effect,  an  appropriate  signal  can  be  used  to  sense  when  it  is  there,  such  as 
maximum  lift  force.  The  UUV  will  then  be  able  to  modify  the  St,  amax,  or  foil  roll 
amplitude  accordingly  to  mimic  the  behavior  from  nature,  if  it  leads  to  better 
efficiency.  The  change  in  signal  strength  will  have  to  be  able  perceived  in  future 
UUVs,  and  will  be  used  in  a  control  feedback  loop  to  keep  the  vehicle  in  ground 
effect,  as  desired. 
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6.  Summary  and  Conclusions 

It  has  been  proven  that  there  is  a  quantifiable  benefit  in  terms  of  lift  and 
thrust  to  be  seen  when  flapping  a  foil  close  to  the  ground,  under  certain 
kinematic  conditions.  Much  work  remains  to  be  done  to  investigate  the  benefit 
and  to  take  advantage  of  it  on  vehicles  that  use  flapping  foils  for  propulsion.  This 
investigation  provided  the  best  possible  starting  point,  as  a  three  dimensional  foil 
was  used  on  an  actuation  system  that  closely  resembles  a  real  underwater 
vehicle. 
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Appendices 


Appendix  1.  MATLAB  code  to  display  results,  "ResultsThesis.m". 
Authored  by  P.  Chierico 

%  Paulstephen  Chierico,  24  OCT  2014 
clc 

%  To  process  the  pitch  bias  sweep  from  "PB240CT" 
pb= [-10:2:10] ; 
pbplot= [-13:2:7] ; 

CLnoGE=-[ 0.8967, 0.8399, 0.7307, 0.5542, 0.4582, 0.3387, 0.1570, 0.017 6, -0.1 127,  .  .  . 
-0.2662,-0.4229] ; 

CLGE=-[ 0.94 86, 0.8123, 0.6947, 0.5589, 0.3850, 0.231 8, 0.1013, -0.0853, -0.24 86, . . . 

-0.4294,-0.6252] ; 
figure ( 1 ) 

CLnoGEeqn=0 . 06785*pb-0. 29011; 

CLGEeqn-0. 07 8354 *pb-0. 2131; 

plot (pbplot, CLnoGE, ' bo ' , pbplot, CLGE, ' rs ' , . . . 

pbplot, CLnoGEeqn, ' b-- ' , pbplot, CLGEeqn, 1 r-- ' , ' LineWidth ' , 2 ) 
grid (gca, 'minor ' ) 

legend('Mean  C_L  f reestream ' , ' Mean  C_L  ground  ef feet Location Northwest ' ) 

%title('Mean  Lift  Coefficient  as  a  function  of  pitch  bias') 

xlabel (' Pitch  Bias  (Ao) ' ) 

ylabel('Mean  CL') 

axis ( [-13  7  -1  0.8] ) 


%  Mean  Coefficient  of  Lift  for  no  GE 
CLnoGEp3=- [0.1355,0.0519,0.0119,0.0295] ; 

CLnoGEp4=-[ 0.1087, 0.07 16, 0.0557, 0.001 6] ; 

CLnoGEp5=-[ 0.0657, 0.1323, 0.1443, 0.0234] ; 

CLnoGEp6=-[ -0.0873, 0.0697, 0.142 6, 0.0858] ; 

%  Mean  Co  of  L  for  GE 

CLGEp3=- [-0.0583,-0.0828,-0.0862,-0.0822] ; 

CLGEp4=- [-0.0048,-0.0193,-0.0573,-0.0869] ; 

CLGEp5=- [-0.0537,0.0368,0.0114, -0.1073] ; 

CLGEp6=-[ -0.12 12, -0.0431, 0.02 64, -0.02 65] ; 

delCLp3=CLGEp3-CLnoGEp3; 
delCLp4=CLGEp4-CLnoGEp4 ; 
delCLp5=CLGEp5-CLnoGEp5; 
delCLp6=CLGEp6-CLnoGEp6 ; 

maxalpha= [20,25,30,35] ; 

figure  (2 ) 

f ig2=plot (maxalpha, delCLp3 ,  ' A  — ' , maxalpha, delCLp4 ,  ' s —  '  ,  .  .  . 

maxalpha, delCLp5 , ' p-- ' , maxalpha, delCLp6 , ' h — ' , ' LineWidth ' ,2) ; 
set (fig2 , { 'markers' } , { 7 ; 5 ; 10 ; 10 } ) 

legend ( ' St=0 . 3 ' , ' St=0 . 4 ' , ' St=0 . 5 ' , ' St=0 . 6 ' , ' Location ' , ' Best ' ) 
grid  on 

%title (' Change  in  Lift  Coefficient  as  a  function  of  Maximum  Angle  of  Attack') 
xlabel  (' \alpha_{max]  (Ao)') 
ylabel ( ' C_L  GE  -  C_L  Free ' ) 
axis (  [19  36  0  0.2] ) 


%  Max  Lift  Force  free 

maxLnoGEp3=- [-2 .9866, -4.2033, -6.4628, -8.2616]; 
maxLnoGEp4=- [-5.2855,-5.3902,-6.2996,-9.3262]; 
maxLnoGEp5=- [-8.4011, -7.8386, -8.2240, -11.2873] ; 
maxLnoGEp6=- [-13.0121,-10.4248, -11.4329,-13. 3272] ; 

%  Max  L  GE 

maxLGEp3=- [-4.7309,-4.8971,-6.6105,-8.6819] ; 
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maxLGEp4=- [-6.3666,-6.1882,-7.4308,-9.8076] ; 
maxLGEp5=- [-9.1873,-8.6612,-9.3484,-12.099] ; 
maxLGEp6=- [-14 .3453, -11 . 6480,-12.1840,-14.3418]; 

delmaxLp3=maxLGEp3-maxLnoGEp3 ;  %  diff 
delmaxLp4=maxLGEp4-maxLnoGEp4 ; 
delmaxLp5=maxLGEp5-maxLnoGEp5 ; 
delmaxLp6=maxLGEp6-maxLnoGEp6 ; 

figure (3) 

f ig3=plot (maxalpha, delmaxLp3 , maxalpha, delmaxLp4 , ' s-- ' , . . . 

maxalpha, delmaxLp5 , 1 p-- ' , maxalpha, delmaxLp6 , ' h-- ' , ' LineWidth ' ,2) ; 
set (fig3, { 'markers' } , { 7 ; 5; 10; 10 } ) 

legend ( ' St=0 . 3 ' , ' St=0 . 4 ' , ' St=0 . 5 ' , ' St=0 . 6 ' , ' Location ' , ' Best ' ) 
grid  on 

%title (' Change  in  Maximum  Lift  Force  as  a  function  of  Maximum  Angle  of  Attack') 
xlabel  (' \alpha_{max]  (Ao)') 

ylabelf'Max  Lift  Force  GE  -  Max  Lift  Force  Free  (N) ') 
axis ( [19  36  0  2] ) 

%  Min  Lift  Force  no  GE 

minLnoGEp3=- [4.9366,5.2389,6.4875,8.7258] ; 
minLnoGEp4=- [6.7889,7.0255,7.7401,9.0427] ; 
minLnoGEp5=- [9.5034, 10.5610, 11.3006, 11.9134] ; 
minLnoGEp6=- [11 . 6982,12.0748,13.7503,15.2993] ; 

%  Min  Lift  Force  GE 

minLGEp3=- [3.5346,4.5951,6.0504,8.1283] ; 
minLGEp4=- [6.0164, 6.5412,7.1570, 9.0405] ; 
minLGEp5=- [8.1668, 9.7931, 10.7058, 11.7424] ; 
minLGEp6=- [11 .4876, 11 . 1195, 13.0805, 14 . 658] ; 

delminLp3=minLGEp3-minLnoGEp3 ; 
delminLp4=minLGEp4-minLnoGEp4 ; 
delminLp5=minLGEp5-minLnoGEp5 ; 
delminLp6=minLGEp6-minLnoGEp6 ; 

figure (4) 

f ig4=plot (maxalpha, delminLp3 , ' ' -- ' , maxalpha, delminLp4 , 1 s — ' , . . . 

maxalpha, delminLp5 , ' p — ' , maxalpha, delminLp6 , ' h — ' , ' LineWidth '  ,  2 )  ; 
set (fig4 , { 'markers' } , { 7 ; 5; 10; 10 } ) 

legend ( ' St=0 . 3 ' , ' St=0 . 4 ' , ' St=0 . 5 ' , ' St=0 . 6 ' , ' Location ' , ' Best ' ) 
grid  on 

%title (' Change  in  Minimum  Lift  Force  as  a  function  of  Maximum  Angle  of  Attack') 
xlabel (' \alpha_{max]  C'o)  ’  ) 

ylabel('Min  Lift  Force  GE  -  Min  Lift  Force  Free  (N) ') 
axis ( [19  36  0  1.5]) 


CTnoGE=-[ -0.2431, -0.317 8, -0.31 96, -0.3005, -0.3142, -0.412, -0.374 8, -0.3775, . . . 
-0.3254,-0.2745,-0.2749]; 

CTGE=- [-0.2483, -0.3255, -0.3446, -0.362,-0.3218,-0.3272,-0.3795,-0.3759, . . . 

-0.3317,-0.2875,-0.2371] ; 
figure (5) 

plot  (pbplot,  CTnoGE,  '  bo--  '  ,  pbplot,  CTGE,  '  rs—  '  ,  '  LineWidth '  ,  2 ) 

%title('Mean  Thrust  Coefficient  as  a  function  of  pitch  bias') 

grid (gca, 'minor ' ) 

xlabel (' Pitch  Bias  ( ^o ) ' ) 

ylabel ( ' C_T ' ) 

legendf'Mean  C_T  f reestream ' , ' Mean  C_T  ground  effect') 
axis ( [-13  7  00.5]) 


%  Mean  Coefficient  of  Thrust  for  no  GE 
CTnoGEp3=- [-0.0978,-0.1030,-0.1258,-0.1131] ; 
CTnoGEp4=- [-0.1128,-0.1868,-0.2120,-0.2326] ; 
CTnoGEp5=- [-0.1327,-0.2475,-0.3403,-0.4016] ; 
CTnoGEp6=- [-0.1479,-0.3697,-0.5100,-0.5777] ; 


%  Mean  Co  of  T  for  GE 

CTGEp3=- [-0.0763, -0.0840, -0.1126, -0.1119] ; 
CTGEp4=- [-0.0947,-0.1878,-0.2263,-0.2153] ; 
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CTGEp5=-[ -0.12 99, -0.24 12, -0.3689, -0.4 190] ; 

CTGEp6=- [-0.1223,-0.3433,-0.5518,-0.5723] ; 

delCTp3=CTGEp3-CTnoGEp3 ; 
delCTp4=CTGEp4-CTnoGEp4 ; 
delCTp5=CTGEp5-CTnoGEp5 ; 
delCTp6=CTGEp6-CTnoGEp6 ; 

figure (6) 

f  ig6=plot  (maxalpha,  delCTp3 ,  1 A--  '  ,  maxalpha,  delCTp4 ,  '  s-  -  1  ,  .  .  . 

maxalpha, delCTp5 , ' p-- ' , maxalpha, delCTp6 , ' h-- ' , ' LineWidth ' , 2 ) ; 
set (f ig6, { 'markers' } , { 7 ; 5; 10; 10 } ) 

legend ( ' St=0 . 3 ' , ' St=0 . 4 ' , ' St=0 . 5 ' , ' St=0 . 6 ' , ' Location ' , ' Best ' ) 
grid  on 

%title (' Change  in  Thrust  Coefficient  as  a  function  of  Maximum  Angle  of  Attack') 
xlabel  (' \alpha_(max]  (Ao) ' ) 
ylabel ( ' CJT  GE  -  C_T  Free ' ) 
axis ( [19  36  -0.03  0.05] ) 

%  Max  Thrust  Force  free 

maxTnoGEp3=- [-2 .217  8, -2 .3124, -2 .6190,-2.7647]; 
maxTnoGEp4=- [-3.1119,-3.4356,-4.0593,-4.3017] ; 
maxTnoGEp5=- [-5.4255,-5.9517,-6.4193,-6.9927] ; 
maxTnoGEp6=- [-8.9669,-8.7929,-9.7747,-10.3269] ; 

%  Max  thrust  force  GE 

maxTGEp3=- [-1 .7600, -2 . 1112, -2 . 3105,-2 .4823]; 
maxTGEp4=- [-2 . 7500,-3.2243,-3.8768,-4.2156] ; 
maxTGEp5=- [-4.7371,-5.3553,-6.2268,-6.4947] ; 
maxTGEp6=- [-8 .2244, -8 .2547, -9. 1805, -9. 8320] ; 

delmaxTp3=maxTGEp3-maxTnoGEp3 ; 
delmaxTp4=maxTGEp4-maxTnoGEp4 ; 
delmaxTp5=maxTGEp5-maxTnoGEp5 ; 
delmaxTp6=maxTGEp6-maxTnoGEp6 ; 

figure  (7 ) 

f ig7=plot (maxalpha, delmaxTp3 , '  A-  -  '  , maxalpha, delmaxTp4 , ' s — ' , . . . 

maxalpha, delmaxTp5 , ' p-- ' , maxalpha, delmaxTp6 , ' h — ' , ' LineWidth '  ,  2 )  ; 
set (fig7 , { 'markers' } , { 7 ; 5; 10; 10} ) 

legend ( ' St=0 . 3 ' , ' St=0 . 4 ' , ' St=0 . 5 ' , ' St=0 . 6 ' , ' Location ' , ' Best ' ) 
grid  on 

%title (' Change  in  Maximum  Thrust  Force  as  a  function  of  Maximum  Angle  of  Attack') 
xlabel (' \alpha_{max]  (Ao)') 

ylabel ('Max  Thrust  Force  GE  -  Max  Thrust  Force  Free  (N) ') 
axis ( [19  36  -0.8  0] ) 

%  Min  Thrust  force  free 

minTnoGEp3=- [0.4861,0.4645,0.3701,0.4076] ; 
minTnoGEp4=-[ 1.04 92, 0.6237, 0.4057, 0.4 163] ; 
minTnoGEp5=- [1 . 6996,1.1352,0.8612,0.6536] ; 
minTnoGEp6=- [4 . 1461,1.9887,1.5649,0.9520]; 

%  Min  Thrust  GE 

minTGEp3=- [0.6656, 0.5932, 0.5787, 0.4189] ; 
minTGEp4=-[ 1.0600, 0.6034, 0.5055, 0.54 91] ; 
minTGEp5=-[ 1.6039, 1.094 6, 0.6837, 0.638 9] ; 
minTGEp6=- [3. 8285, 1 .8215,1.4235,1.2111]; 

delminTp3=minTGEp3-minTnoGEp3 ; 
delminTp4=minTGEp4-minTnoGEp4 ; 
delminTp5=minTGEp5-minTnoGEp5 ; 
delminTp6=minTGEp6-minTnoGEp6 ; 

figure (8) 

f ig8=plot (maxalpha, delminTp3 maxalpha, delminTp4 , ' s — ' , . . . 

maxalpha, delminTp5 , ' p-- ' , maxalpha, delminTp6 , ' h-- ' , ' LineWidth ',2); 
set (f ig8, ( 'markers' } , { 7 ; 5; 10; 10 } ) 

%title (' Change  in  Minimum  Thrust  Force  as  a  function  of  Maximum  Angle  of  Attack') 
legend ( ' St=0 . 3 ' , ' St=0 . 4 ' , ' St=0 . 5 ' , ' St=0 . 6 ' , ' Location ' , ' Best ' ) 
grid  on 

xlabel (' \alpha_[max]  (Ao) ' ) 

ylabel ('Min  Thrust  Force  GE  -  Min  Thrust  Force  Free  (N) ') 
axis ( [19  36  -0.3  0.4] ) 
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St= [0.3:0. 1:0. 6] ; 

delCL30= [delCLp3 (1,3) ,delCLp4 (1,3) , delCLp5 ( 1 , 3) , delCLp6 ( 1 , 3) ] ; 
delminL30= [delminLp3 (1,3) , delminLp4 (1,3) , delminLp5 (1,3) , delminLp6 (1,3) ] ; 
delmaxL30= [delmaxLp3 (1,3) , delmaxLp4 (1,3) , delmaxLp5 (1,3) , delmaxLp6 (1,3) ] ; 
delCT30=[delCTp3 (1,3) ,delCTp4 (1,3) , delCTp5 ( 1 , 3) , delCTp6 ( 1 , 3) ] ; 
delminT30= [delminTp3 (1,3) , delminTp4 (1,3) , delminTp5 (1,3) , delminTp6 (1,3) ] ; 
delmaxT30= [delmaxTp3 (1,3) , delmaxTp4 (1,3) , delmaxTp5 (1,3) , delmaxTp6 (1,3) ] ; 

figure (9) 

f ig9=plot  (St,  delCL30,  '  bv--  '  ,  St,  delCT30,  '  rx—  '  ,  '  LineWidth '  ,  2 )  ; 
set (fig9, { 'markers ' } , {7; 10} ) 

%title (' Change  in  Mean  C_L  or  Mean  C_T  as  a  function  of  St  for  \alpha_{max}=30Ao ' 
grid  on 
xlabel  ( 'St' ) 

ylabel ( ' \Delta  Mean  C_L  or  \Delta  Mean  C_T ' ) 

legend (' \Delta  Mean  C_L','\Delta  Mean  C_T ',' Location ',' Best ' ) 
axis  (  [0.28  0.62  -0.02  0.14]) 

figure (10) 

plot (St, delminL30 , ' gv-- ' , St , delmaxL30 , ' mA-- ' , ' LineWidth ' , 2 ) 

%title (' Change  in  Min  L  or  Max  L  as  a  function  of  St  for  \alpha_{max}=30Ao ' ) 
grid  on 
xlabel  (  'St' ) 

ylabel ( ' \Delta  Min  L  or  \Delta  Max  L  (N) ' ) 

legend (' \Delta  Min  L','\Delta  Max  L ',' Location ',' Best ' ) 

axis ( [0.28  0.62  0  1.4]) 

figure ( 11 ) 

plot (St, delminT30 , ' gv-- ' , St , delmaxT30 , ' mA-- ' , ' LineWidth ' , 2 ) 

%title (' Change  in  Min  T  or  Max  T  as  a  function  of  St  for  \alpha_{max}=30 Ao ' ) 
grid  on 
xlabel ( 'St' ) 

ylabel (' \Delta  Min  T  or  \Delta  Max  T  (N) ' ) 

legend (' \Delta  Min  T','\Delta  Max  T ',' Location ',' Best ' ) 

axis ( [0.28  0.62  -0.6  0.2]) 

%  Section  for  contour  plots 


a=0 . 3 : 0 . 1 : 0 . 6; 
b=20 : 5 : 35; 

[X, Y]  =  meshgrid (a, b) ; 

ZCLnoGE= [CLnoGEp3; CLnoGEp4 ; CLnoGEp5; CLnoGEp6] ' ; 
figure (12) 

contour (X, Y, ZCLnoGE, ' ShowText ' , ' on ' , ' LineWidth ' , 3) 

%title('Mean  Lift  Coefficient  in  Freestream') 

xlabel ( ' St ' , ' fonts ize ' , 12 ) 

ylabel ( ' \alpha_{max}  ( Ao) ' , ' fontsize ' , 12) 

grid  on 

axis  square 

ZCLGE= [CLGEp3 ; CLGEp4 ; CLGEp5 ; CLGEp6 ] ' ; 
figure (13) 

contour (X, Y, ZCLGE, ' ShowText ' , ' on ' , ' LineWidth ' , 3) 

%title('Mean  Lift  Coefficient  in  Ground  Effect') 

xlabel ( ' St ' , ' fontsize ' , 12 ) 

ylabel ( ' \alpha_{max}  ( Ao) ' , ' fontsize ' , 12) 

grid  on 

axis  square 

ZdCL=ZCLGE-ZCLnoGE; 
figure (14) 

contour (X, Y, ZdCL, ' ShowText ' , ' on ' , ' LineWidth ' , 3) 

%title (' Difference  in  Mean  Lift  Coefficient  (GE  -  Free)') 

xlabel ( ' St ' , ' fontsize ' , 12 ) 

ylabel ( ' \alpha_{max}  ( Ao) ' , ' fontsize ' , 12) 

grid  on 

axis  square 

ZCTnoGE= [CTnoGEp3 ; CTnoGEp4 ; CTnoGEp5 ; CTnoGEp6 ] ' ; 
figure (15) 

contour (X, Y, ZCTnoGE, ' ShowText ' , ' on ' , ' LineWidth ' , 3) 
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%title('Mean  Thrust  Coefficient  in  Freestream') 

xlabel ( ' St ' , ' fonts ize ' ,12) 

ylabel ( ' \alpha_{max}  ( Ao) ' , ' fontsize ' ,12) 

grid  on 

axis  square 

ZCTGE= [CTGEp3 ; CTGEp4 ; CTGEp5 ; CTGEp6 ] ' ; 
figure (16) 

contour (X, Y, ZCTGE, ' ShowText '  , ' on ' , ' LineWidth ' , 3) 

%title('Mean  Thrust  Coefficient  in  Ground  Effect') 

xlabel ( ' St ' , ' fontsize 1 ,12) 

ylabel ( ' \alpha_{max}  ( Ao) ' , ' fontsize ' ,12) 

grid  on 

axis  square 

ZdCT=ZCTGE-ZCTnoGE; 
figure (17) 

contour (X,Y, ZdCT, 'ShowText' , 'on', 'LineWidth' ,3) 

%title (' Difference  in  Mean  Thrust  Coefficient  (GE  -  Free)') 

xlabel ( ' St ' , ' fontsize ' , 12 ) 

ylabel ( ' \alpha_{max}  ( Ao) ' , ' fontsize ' , 12) 

grid  on 

axis  square 
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Appendix  2.  MATLAB  code  for  phase  averaging,  "PhaseAvgO.m". 
Authored  by  Professor  S.  Licht 


function  [phase_time,  phase_data]  =  PhaseAvgO (  timep,  pdata,  time,  data,  frequency  ) 
%Phase  averages  FMdata  over  one  2pi  cycle 

%Find  cycle  start  using  peak  roll  position: 
cycle_to_start  =  2;  %discard  first  peak 

%Use  roll  position  peaks  to  estimate  start  of  cycles 

[ppy,  ppx]  =  findpeaks (pdata) ; 

time_step  =  mean (diff (time) ) ; 

period  =  mean (diff (timep (ppx) )) ; 

roll_mean  =  mean (pdata); 

roll_amplitude  =  (max (pdata) -min (pdata) ) /2 ; 

%Use  evenly  spaced  values  at  mean  sampling  frequency: 
samples_per_cycle  =  period/time_step; 

phase_time  =  linspace (O,period-time_step, samples_per_cycle) ; 

for  i  =  1 : (length (ppx) -cycle_to_start) 

%Use  the  peaks  in  the  position  data  to  identify  approximate  start  and 
%end  of  this  cycle  in  the  position  data: 
start_index  =  ppx (i+cycle_to_start-l ) ; 
end_index  =  ppx (i+cycle_to_start) ; 

time_past_peak  =  timep (start_index:end_index)  -  timep (start_index) ; 
position  =  pdata (start_index:end_index) ; 

%Fit  a  sinusoid  to  the  data: 

fit  =  @(b,x)  b  (1)  .*  (cos  (2*pi*x. /period  +  b(2)))  +  b(3);  %  Function  to  fit 

fen  =  @ (b)  sum ( (f it (b, time_past_peak)  -  position)  .  ~2 )  ;  %  Least- 

Squares  cost  function 

s  =  fminsearch (fen,  [roll_amplitude;  0;  roll_mean] ) ;  % 

Minimise  Least-Squares  with  starting  guesses 

time_shift  =  period  *  s(2)/(2*pi);  %this  is  the  only  output  used  from  the  fit 

%Now  use  the  time  shift  estimate  to  identify  the  start  and  end  of  this 
%cycle  in  the  data  that  needs  to  be  phased  averaged: 
cycle_start_time  =  timep (ppx (i+cycle_to_start-l )) -time_shift; 
if  (min (size (data) )  ==  1) 

resized_data (i, :)  =  interpl (time-cycle_start_time, data, phase_time) ; 

else 

for  k  =  1 : size (data, 2 ) 

resized_data (i, : , k)  =  interpl (time-cycle_start_time, data (:, k) , phase_time) 

end 

end 

end 

%Test  Comparison 
%figure (101) 

%plot (phase_time, resized_data) 

if  (min (size (data) )  ==  1) 

averaged_data  =  mean (resized_data) ; 

else 

for  k  =  1 : size (data, 2 ) 

averaged_data (k, : )  =  mean (squeeze (resized_data (:,:, k) ) )  ; 

end 

end 

phase_data  =  averaged_data; 
end 
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Appendix  3.  MATLAB  code  for  finding  force  and  moment  data,  "Parse2FM0.m". 


Authored  by  S.  Rauworth,  modified  by  P.  Chierico 

function  [  FMdata,  TVPdata2,  Fldata,  F2data  ]  =  Parse2FM0 (  FSdata,  TVPdata,  num  ) 
%Takes  parsed  force  sensor  data,  user  chooses  good  data  from  plot, 

%Force/Moment  data  outputted  in  x,  y,  and  z  directions 

%  new  ca  matrices  due  to  results  from  hanging  directly  off  bearings 
CmSl= [484.6,  8.826,  -28.18;  9.719,  432.8,  35.28;  39.23,  -72.35,  1207]; 

CmS2= [481.8,  -33.13,  -64.44;  -6.871,  494.0,  -54.16;  29.54,  -70.20,  1109]; 

time  =  FSdata ( : , 1 ) ; 

plot  (FSdata (:, 6) ) 

[x,y]  =  ginput  (2)  ; 

data  (1,1)  =  round  (x(l)); 

data  (1,2)  =  round(x(2)); 

%  2  =  xl  and  4  =  yl  and  6  =  zl 

S1V  =  [FSdata (:, 2)  FSdata(:,4)  FSdata(:,6)  ]'; 

S1F  =  CmSl*SlV; 

%  3  =  x2  and  5  =  y2  and  7  =  z2 

%  new  change,  leave  everything  in  sensor  frame,  until  forces  are  added 
S2V  =  [FSdata (:, 3)  FSdata (;, 5)  FSdata (:, 7 )]' ; 

S2F  =  CmS2*S2V; 

m  =  data ( 1 ) ; 
n  =  data  (2 )  ; 

for  i  =  m:m+10 

if  mod (i, 10)  ==  0 
g  =  i; 

end 

end 

fl  =  (g/10) ; 
bl  =  (f 1*10) -9; 

for  j  =  n:n+10 

if  mod ( j , 10)  ==  0 

h  =  j; 

end 

end 

f2  =  (h/10)  ; 
b2  =  ( f 2  *10) —  9 ; 

TVPdatanew  =  TVPdata (f 1 : f2 , 2 : 7 ) ; 

TVPtimenew  =  0  :.  05  :  (length (TVPdatanew) /20) -. 05; 

TVPdata2  =  [TVPtimenew',  TVPdatanew] ; 

FS1  =  S1F ( : , bl : b2 ) ; 

FS2  =  S2F  ( : ,bl :b2) ; 

FS  =  FS1+FS2;  %  what  is  this  doing? 

time2  =  0 :. 005 :(( length (FS1 )) /200 )-. 005 ; 

FMdata  =  [time2 ' ,  FS ' ] ; 

Fldata  =  [time2 ' ,  FS1 ' ] ; 

F2data  =  [time2 ' ,  FS2 ' ] ; 

end 
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Appendix  4.  MATLAB  code  for  filtering,  "BFiltO.m". 


Authored  by  S.  Rauworth,  modified  by  P.  Chierico 

function  [  xout,  yout,  zout  ]  =  BFiltO (  xdata,  ydata,  zdata) 

%Moving  Average  Filter 

%Filters  force  sensor  data  forwards  and  backwards  with  moving  average. 
[b,a]  =  butter  (5, 0 . 1 low ') ; 

xl  =  xdata ; 
x2  =  ydata ; 
x3  =  zdata ; 

xout  =  filtfilt (b, a,xl) ; 
yout  =  f iltf ilt (b, a, x2 ) ; 
zout  =  f iltf ilt (b, a, x3) ; 

end 
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Appendix  5.  MATLAB  code  for  analyzing  the  data,  "AnalysisO.m". 


Authored  by  S.  Rauworth,  modified  by  P.  Chierico 

%  Paulstephen  Chierico,  modified  code  originally  by  Sam  Rauworth 

%Use  this  script  to  fully  parse,  organize,  and  analyze  data  from  individual 

%tests 

close  all 
clear  all 
clc 

addpath ( ' PB240CT ' ) 
addpath ( ' GEBias240CT ' ) 

f ilenm= ' TEST1 6L . lvm ' ; 
testnum=146;  %  used  up  to  146 
Fr=l . 65; 

AMX= 12; 

AMY=27; 

St=0 . 6 ; 
maxaoa=35; 

[ TVP ,  FSd]  =  DataParseO (f ilenm) ;  %This  function  parses  the  data  from  the  above  file 
%TVP  is  Torque  Velocity  Position  data  from  Control  Card 
%FSd  is  Force  Sensor  data 

f  ileB=  '  ZB1 6L  .  lvm '  ; 
f ileE= ' ZE1 6L . lvm ' ;  %  R4GE 

LDMf ile= ' LDM1 6L . lvm ' ; 

Vavg  =  VfunO (LDMf ile, testnum) ;  %Change  the  file  name  here  for  the  laser  range  data  to 
find  average  velocity  of  the  test. 

%User  is  asked  to  choose  range  of  velocity  to  average 

[FSdz]  =  zeroingO (FSd, fileB, fileE) ;  %Finds  the  Force  Sensor  data  after  zeroing  it. 

[FMd,  TVPd,  Fid,  F2d]  =  Parse2FM0 (FSdz,  TVP,  testnum);  %Finds  the  Force  and  Moments 
outputted  by  the  force  sensors. 

%User  is  asked  to  choose  the  data  range  which  is  output  to  Test_Range.txt 
%FMd  is  the  Force  Moment  data 

%TVPd  is  the  TVP  data  corresponding  to  the  chosen  FMd 

FMdx  =  FMd(:,l);  %Sets  variables  for  individual  force  directions 

FMdy  =  FMd ( : , 2 ) ; 

FMdz  =  FMd ( : , 3 ) ; 

[x,y,z]  =  BFiltO ( FMdx, FMdy , FMdz ) ;  %Applies  butterworth  filter  to  those  forces 

%  To  switch  from  sensor  frame  forces  into  the  body  frame  of  the  fin 
Fldbody= [-Fid ( : , 2) , -Fid ( : , 3) , -Fid ( : , 4) ] ;  %  taking  time  out  -  lv  first  column 
F2dbody=[F2d(:  ,2)  ,F2d(  :  ,3)  ,-F2d( :  ,4)  ]  ;  % 

%  torque  ABOUT  the  axis 
%  so  xt  only  contribution  is  y-force 
%  yt  only  contributions  are  x-force  and  z-force 
%  zt  only  contribution  is  y-force 

xtorque  =  Fldbody ( : , 2 ) * (tO . 04 61 ) tF2dbody ( : , 2 ) * (tO . 04 61 ) ; 

ytorque  =  Fldbody (:, 1 ) * (-0 . 04 61 ) tF2dbody ( : , 1) * (-0 . 04 61 ) tFldbody ( : , 3 ) * (- 

0.0349) tF2dbody( : ,3) * (0.0572) ; 

ztorque  =  Fldbody ( : , 2 ) * (0 . 034 9 ) +F2dbody ( : , 2 ) * ( -0 . 0572 ) ; 

[Tx,Ty,Tz]  =  BFiltO (xtorque, ytorque, ztorque) ; %Applies  butterworth  filter  to  those 
torques 

FS_time  =  (FMd(:,l));  %Time  used  with  forces  and  moments  (different  frequency  than 
control  card) 

TVP_time  =  (TVPd(:,l));  %Time  used  with  TVP  (different  frequency  than  force  sensors) 

%  %  Position  (sanity  check) : 

%  [phase_time.  Phased?]  =  PhaseAvgO (TVP_time,  TVPd(:,5),  TVP_time,  TVPd(:,5),  Fr) ; 
%Can  change  variable  name  and  2nd  function  input  to  phase  average  different  variables 
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%  plot (phase_time, PhasedP) 

%  title (' Position  vs.  Phase  Sanity  Check') 

[phase_time,  PhasedTy]  =  PhaseAvgO (TVP_time,  TVPd(:,5),  FS_time,  Ty,  Fr)  ; 
[phase_time,  PhasedTz]  =  PhaseAvgO (TVP_time,  TVPd(:,5),  FS_time,  Tz,  Fr)  ; 

Lift  =  PhasedTy/O . 435;  %Converts  phase  averaged  torques  into  forces  by  dividing  out 

moment  arm  or  distance  from  axis  of  rotation  to  center  of  pressure 

Thrust  =  PhasedTz/O . 435 ; %Min  CoP  is  .155m,  max  is  .555m,  70  percent  is  .435m 


for  i  =  1 : length (Thrust)  %  The  below  removes  the  NaN  values  from  the  thrust  and  lift 
matrices 

if  isnan (Thrust (i) )  ==  0 

Thrust2 (i, 1) =Thrust  (i) ; 

end 

end 

for  i  =  1 : length  (Lift) 

if  isnan (Lift  (i) )  ==  0 
Lift2 (i, 1) =Lift  (i) ; 

end 

end 

%  To  find  the  mean  thrust  coefficient 

Ct  =  (2  *mean (Thrust2 ) ) / (1000* (VavgA2 ) * . 1* . 4 ) ; 

Cl  =  (2  *mean (Lift2) ) / (1000* (Vavg^2 ) *  .  1  *  .  4 )  ; 

MinLift=min (Lift) ; 

MaxLift=max (Lift) ; 

MinThrust=min (Thrust) ; 

MaxThrust=max (Thrust) ; 


%  Change  file  name  and  structure  name  for  different  testing  sets 
load ( ' PstTests .mat ' ) 

Struct (testnum) . name  =  filenm; 

Struct (testnum) . St  =  St; 

Struct (testnum) .maxaoa  =  maxaoa; 

Struct (testnum) . Freq  =  Fr; 

Struct (testnum) .AMX  =  AMX; 

Struct (testnum) .AMY  =  AMY; 

Struct (testnum) .Raw_FS  =  FSd; 

Struct (testnum) .Raw_TVP  =  TVP; 

Struct (testnum) . ChosenFd  =  FMd; 

Struct (testnum) . ChosenFl  =  Fid; 

Struct (testnum) . ChosenF2  =  F2d; 

Struct (testnum) . ChosenFY  =  FMdy; 

Struct (testnum) .ChosenFX  =  FMdz; 

Struct (testnum) . ChosenTVP  =  TVPd; 

Struct (testnum) . TorqueY  =  ytorque; 

Struct (testnum) . TorqueZ  =  ztorque; 

Struct (testnum) . FilteredFY  =  y; 

Struct (testnum) . FilteredFZ  =  z; 

Struct (testnum) . FilteredYtorque  =  Ty; 

Struct (testnum) . FilteredZtorque  =  Tz; 

Struct (testnum) . PhaseAvgY  =  PhasedTy; 

Struct (testnum) . PhaseAvgZ  =  PhasedTz; 

Struct (testnum) . Lift  =  Lift; 

Struct (testnum) . Thrust  =  Thrust; 

Struct (testnum) . ForceTime  =  FS_time;  %  time 
Struct (testnum) . PollTime  =  TVP_time;  %  time2 
Struct  (testnum)  . ThrustCoeff  =  Ctj- 
Struct  (testnum)  .  LiftCoeff  =  Cl; 

Struct (testnum) .MinL  =  MinLift; 

Struct (testnum) .MaxL  =  MaxLift; 

Struct (testnum) .MinT  =  MinThrust; 

Struct (testnum) .MaxT  =  MaxThrust; 

save  ( ' PstTests .mat 1 , ' Struct ' , 1 -append' ) 
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